MATHEMATICAL STUDY OF ROTATING FLUIDS 
WITH RESONANT SURFACE STRESS 
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ANNE-LAURE DALIBARD AND LAURE SAINT-RAYMOND 



Abstract. We are interested here in describing the linear response 
of a highly rotating fluid to some surface stress tensor, which admits 
fast time oscillations and may be resonant with the Coriolis force. In 
addition to the usual Ekman layer, we exhibit another - much larger 
- boundary layer, and we prove that for large times, the effect of the 
surface stress may no longer be localized in the vicinity of the surface. 
From a mathematical point of view, the main novelty here is to introduce 
some systematic approach for the study of boundary effects. 



The goal of this paper is to understand the influence of a surface stress - 
depending on time - on the evolution of an incompressible and homogeneous 
rotating fluid. More precisely, we are interested in the effects of a resonant 
forcing, i.e. of a stress oscillating with the same period as the rotation of 
the fluid. 

In the non-resonant case, the works by Desjardins and Grenier [5] then by 
Masmoudi [16] show that the wind forcing creates essentially some boundary 
layer in the vicinity of the surface, which contributes to the mean motion 
by a source term, known as the Ekman pumping. For a precise description 
of the method leading to such convergence results, we refer to the book [3] 
by Chemin, Desjardins, Gallagher and Grenier. 

Here the situation is much more complicated since the resonant part of the 
forcing will be proved to generate another boundary layer with a different 
typical size, and may overall destabilize the whole fluid with the apparition 
of a vertical profile. We give here a precise description of these (linear) 
effects of the Coriolis force in presence of resonant wind. 



Let us first present the mathematical framework of our study. 

1.1. A linear model for rotating fluids. 

• Our starting point is the linear version of the homogeneous incom- 
pressible Navier-Stokes system in a rotating frame 



1. Introduction 




d t u + Vp = T + u A n 
V • u = 
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where T denotes the frictional force acting on the fluid, O is the rotation 
vector, and p is the pressure defined as the Lagrange multiplier associated 
with the incompressibility constraint. We assume that equation is al- 
ready in a nondimensional form, meaning that all unknowns and parameters 
are dimensionless. For a precise dimensional analysis, we refer for instance 
to [H] (section 1.3). 

We assume further that the rotation vector 0, is constant, homogeneous, 
and has constant vertical direction, which we denote by e^. Moreover, we 
wish to study the limit of fast rotation, i.e. — ► oo. Hence, we set 

$7 := — e3, with e — > 0, 
e 

where the parameter e is called the Rossby number. 

• We consider the motion in some horizontal strip 

uj = u h x [0, 1] 

where the bottom and upper surface of the fluid are assumed to be flat at 
z = and z = 1. For the sake of simplicity, we restrict our attention to the 
case when uJh = T 2 is the the two-dimensional torus. 

As boundary conditions on the upper surface, we enforce 

(1.2) = °' 

d z u h \ z=l = (3a e , 

where (5 is a positive constant and a e is a given stress tensor of order one, 
describing the stress on the surface of the fluid. 

At the bottom we use the Dirichlet boundary condition 

(1.3) u lz=0 = 0. 

• At last, we assume that frictional forces T are given by 

T = A h u + ud zz u, 

such a choice is classical in the rotating fluids literature, see for instance 
[H HH [T7]. We refer to paragraph 16.11 for an attempt of justification in a 
geophysical context. 

Hence, our goal is to study the asymptotic behaviour as e — ► of the 
solution of 

d t u + -e 3 Au + Vp- A h u - vd zz u = , 

(1.4) e 

V • u = 0, 



supplemented with the boundary conditions (jl.2p - (jl.3p . depending of the 
order of magnitude of the vertical viscosity v. 
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1.2. Formal study of the asymptotics. 

The system (|1.2j) - (ll.4p has already been studied by several authors, see 
for instance [16^15]. and also [HE] when Dirichlet boundary conditions are 
enforced at the top and at the bottom. Before describing the precise issues 
we wish to study in the present paper, let us recall briefly some of the main 
results and techniques for singular perturbation problems. 

• The first step is to determine the geostrophic motion. The only way 
to control the Coriolis force as e — ► is to balance it with the pressure 
gradient term (see for instance [Mj). Hence in the limit, A u must be a 
gradient 

(1-5) e 3 A iCt an = -Vp 

which leads to 

int _ Y7-L 

""mean v h f 

where the limit pressure and thus the limit velocity are independent of z. 
In particular, tt™ean ^ s a two-dimensional, horizontal, divergence-free vector- 
field. The fluid being limited by rigid boundaries, from above and below, 
the divergence-free condition leads indeed to U3 = (at least to first order 
in e). In other words, all the particles which have the same Xh have the 
same velocity. The particles of fluid move in vertical columns, called Taylor- 
Proudman columns. That is the main effect of rotation and a very strong 
constraint on the fluid motion. 

As the domain evolution is limited by two parallel planes, the height of 
Taylor-Proudman columns is constant as time evolves, which is compatible 
with the incompressibility constraint. We can then prove that the columns 
move freely and in the limit of high rotation the fluid behaves like a two- 
dimensional incompressible fluid. Integrating the motion equation (jl.4j) with 
respect to z and taking formal limits as e — > leads indeed to 

(1.6) . , 

v h ' a mean ~ u • 

Note however that on the boundary of the domain, where the velocity is 
prescribed, the z independence is violated. That leads to vertical bound- 
ary layers modifying the limit equation (jl.6p . which will be investigated in 
the rest of the paper. 

• Before starting with the precise study of these boundary layers, let us 
now describe what happens for the three-dimensional ageostrophic part 
of the initial data, i.e. the part of the initial data that does not satisfy the 
geostrophic constraint (jl.5p . The dominant process is then governed by the 
Coriolis operator 



(1.7) 



L : u e V h-> P(e 3 Au)e V , 
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where Vq denotes the subspace of L 2 (u>) of divergence-free vector fields hav- 
ing zero flux both through the bottom and through the surface 

V = {u G L 2 ([0, 1] x T 2 ) / V • u = and u 3 | 2=0 = u 3 \ z=1 = 0}, 

and P denotes the orthogonal projection onto Vo in L 2 (uj). Notice that in 
general, Vq is strictly smaller than the space of divergence-free vector fields 
in L 2 (lo), and consequently P is different from the Leray projector. 
The equation 

edtu + Lu = 

turns out to describe the propagation of waves, called Poincare waves. More 
precisely, one can prove (see for instance p3] , [I] and Appendix A at the end 
of this paper for more details) that there exists a hilbertian basis of Vo, de- 
noted by (A r fc) fcgZ 3\ v { } ) constituted of eigenvectors of the linear penalization: 
for all k G Z 3 \ {0}, we have 

(1.8) LN k = F(e 3 AN k ) = i\ k N k , where A fc = ^ 

^J\k h \ 2 + (irk 3 ) 2 

That means that the three-dimensional part of the initial data generates 
waves, which propagate very rapidly in the domain (with a speed of order 
e _1 ). The time average of these waves vanish, like their weak limit, but they 
carry a non-zero energy. 

1.3. Resonant forcing. In view of the remarks of the previous paragraph, 
it seems interesting, in order to study possible resonances between the sur- 
face stress and the Coriolis operator L, to consider in (jl.2p a stress tensor 
of the form 

a e (t,x h ) = a f 

with a G L°°([0, oo) x T 2 ) almost periodic in its first variable, i.e. 

JfiT e ik h -x h 

' VP) "-ft; 1 
k h ez 2 Ate A/ 

where M is a finite set. The corresponding boundary layer terms are then 
expected to oscillate with the frequencies [x/e, with either /i£Mor/i = —\ k 
for some k G Z 3 . The construction of such boundary layer terms is relatively 
well understood (see for instance [U [T6| I17j). insofar as [i ^ ±1. When 
\fi | = 1, the classical construction of boundary layers fails; the usual way 
to get round this difficulty is to assume that the initial data and the stress 
tensor satisfy some spectral assumptions, in order to avoid the apparition 
of the frequencies fj, = ±1 altogether. 

Our goal in this paper is precisely to study the influence of such res- 
onant frequencies on the global behaviour of the fluid, starting with the 
boundary layers. To that end, we have developed a systematic way of com- 
puting the boundary layer profiles associated with some given boundary 
conditions; our main result in that regard is stated in the next paragraph, 



(1.9) a (T,x h )= ^ ^£v>AK r "< 
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and proved in section [3j Next, we use the boundary layer profiles so defined 
in order to construct an approximate solution for equation (|1.4p . supple- 
mented with (|1.2j )-l fL3]) . and we prove a strong convergence result for (jl.4[) . 

2. Main results 

2.1. Description of the boundary layers. 

We begin with the construction of boundary layers. Let us first emphasize 
that since equation (|1.4p is linear, we can work with a finite number of 
Fourier modes in the horizontal domain and in time. Note that on the 
contrary, because of the boundary conditions at z = and z = 1, there is a 
strong coupling between the vertical modes. 

Hence, let N > be an arbitrary integer, and let Mq be a finite set such 
that M C Mq. We consider some arbitrary boundary conditions 5® and 5\ 
which take the form 

(2.1) Si(r,x h )= E Si(^kh)e^ T e ik ^\ j = or 1. 

\k h \<Nn&M 

Here and in the whole paper, the superscript (resp. 1) stands for func- 
tions associated with some boundary conditions at the bottom (resp. at the 
surface) . 

Our goal is to construct some stationnary boundary layer profiles, denoted 
by v®,v l , which have respectively exponential decay with respect to z and 
1 — z, are exact solutions of equation (|1.4j) . and satisfy 

v° h{z=0 {t,x h ) = S° h (-x h ), 

(2 - 2) ft \ 

d z v h \ z=Q {t,x h ) = S h [-,Xhj ■ 

Notice that we do not enforce boundary conditions on both sides for vi, and 
that we do not specify the boundary condition on the vertical component 
of each function v 3 : indeed, the vertical component of v 3 is dictated by the 
assumption that v 3 is divergence free and that its dependance on the vertical 
variable z is given by a decaying exponential. Similarly, the trace of v° at 
z = 1 is imposed by the exponential profile condition. At last, we do not 
specify any initial data for v ^ 1 , for the same reasons as above; we only 
require that [|«u_p||.£,2 = o(l) as e,u — > 0. 

However that construction fails if some particular coefficients S 3 (kh,n) in 
the boundary condition are not identically zero (see Remark on page 15). 
This leads to the following definition: 

Definition 2.1. Assume that the boundary conditions Si are given by 

*h) = E E %0h k h )e^e^ . 
k h fj. 
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We define the resonant part S h TP!i of the boundary conditions by 

1 

|^(1,0/ ' 

We will say that a boundary condition 5 3 h is non-resonant if Si = 0. 

In the resonant case, we will indeed see that the boundary profiles are not 
stationnary. More precisely, we will prove the following result 

Theorem 2.2. Let 5^,5^ be given by (|2.ip . Then there exist v ^ 1 which 
are exact solutions of (|1.4|) supplemented with (|2.2j) . and such that v° decays 
exponentially with z, and v 1 with 1 — z. Moreover, each function v 3 (j = 
or 1) can be written as 

v 3 = v 3 + v 3 + v 3 res 

where the stationnary boundary profiles v 3 , v 3 satisfy the following estimates 
(2.3) 



\rf l h°°(R+,L a (u,)) + -F=lK[|L«>(R+,La(w)) < C (^) + ^ PV 
V ev 

1 1+2 j 

sJil 1 / e + V^V I, -j,. ™ x 



J h\h 



KIU«(R+,La(«)) + ^ — j ||w 3 ||l°°(R+,L2H) < C I — py= 

w;/iife i/ie resonant part v 3 re s satisfies 

(2.4) Vt>0, II^WII^H^^^H^JI, < S)3 = 0, 
w/iere 

11*1= E E $(/*.*fc)i a 

M6M |%|<jV 

and C is a nonnegative constant depending on N . 

Theorem 12.21 will be proved in section [3j The definition of the boundary 
layer operator B is then as follows: 

Definition 2.3. Let 5®, 5\ be given by (|2.1|) . We denote by B the bilinear 
operator such that with the notations of Theorem \2.2[ 



Remark 2.4. (i) As we shall see in the course of the proof, the terms v 3 
correspond to the usual Ekman layers, for which the typical size of the 
boundary layer is \JTv. The corresponding boundary conditions are given by 

K(T>*h)= E E K^ k h)^^ kh ' Xh - 

\k K \<N\p\*\ 
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On the contrary, the terms v 3 are due to the quasi-resonant modes, for 

which = \ and kh 7^ 0; for these modes, the typical size of the boundary 
layer is much larger, of order yfev/(yfe + (ez^) 1//4 ). 

(ii) The last terms v 3 res are due to resonant forcing on the modes = 
l,kh = 0. Notice that for these modes, the estimate is not global in time: 
indeed, the typical size of the boundary layer is \J~ui. 



5{{r,x h )= 

In particular, for large times (t 3> v~ x )> the boundary layer penetrates the 
interior of the fluid. 

(iii) As outlined above, the boundary layer term v° (resp. v 1 ) does not 
vanish on z = 1 (resp. on z = 0). Precisely, we find that there exists a 
positive constant C ( depending on N and M ) such that 

C \\ ,1 nr.^( C 



v? z= i = O ( exp ( ^= ) ) , vl =0 = O ( V^exp 



res 1 2: 



SfU = O (exp (-J^)) . »|U = O (W' 4 exp (-^ 



2.2. Construction of approximate solutions to (|1.4|) - (|1.2p(|1.3|) . 

Once the mechanism of construction of boundary layers is understood, 
one possible application lies in the definition of an approximate solution of 
equation (jl.4p . with a view to derive a limit system for this equation. This 
approximate solution is the sum of boundary terms u BL , obtained as above, 
and interior terms u mt . 

• Hence, we now explain the asymptotic behaviour of the interior part 
of the solution. Following the multi-scale analysis initiated in the previous 
paragraph, we expect the solution u e to (jl.4j) to behave like some function 
exp(— tL/e)u™*(t), where L is the Coriolis operator defined by (jl.7)l . 

In order to understand the evolution with respect to the slow time vari- 
able, the idea is then to get rid of the penalization term by filtering out the 
oscillations in equation (II. 4h (see [H1I23]), that is, by composing equation 
(II. 4j) by the Coriolis semi- group exp(tL/e). 

The filtered function u £i i(t) := exp(tL / e)u e (t) satisfies a linear equation 
with vanishing viscosity (and without any penalization term); passing to the 
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where SEkman '■ Vq — * Vo is a linear, positive and continuous operator re- 
sulting from the non commutation between the vertical Laplacian z/A 2 with 
boundary conditions and the Coriolis semi-group (see [1] and (|5.2ip below 
for a precise definition). 

• The approximation of the function u t constructed in this paper is actu- 
ally much more precise than the mere function exp(— tLjt)v^jf . Indeed, we 
will need to build boundary and corrector terms, which are all small in 
L 2 norm, and thus do not play a role in the final convergence result, but are 
necessary in order that equation (II. 4h is approximately satisfied. 

2.3. Convergence result. 

Theorem 2.5. Let 7 G Vo, and let a be given by (|1.9[) . Let u e G C(R + , Vb)n 
L 2 0C (R + , H 1 (T 2 x [0,1])) be the unique solution of (|1.4j) supplemented with 
([rj ]) - ([T3l) . and let uf f G C(R+, V" )nL 2 (R + , H\(T: 2 x [0, 1])) be the solution 
of equation (|2.5p . 

Assume that a has a finite number of Fourier modes, i.e. a satisfies (|2.ip . 
Then under the technical scaling assumption (|4.17p on the parameters e, v 
and (5, we have, as e, v — > 0, 



Remark 2.6. (i) That result extends previous works by Masmoudi |16| and 
Chemin, Desjardins, Gallagher and Grenier [1]. They have indeed studied 
analogous boundary problems for rotating fluids, but have used in a crucial 
way a spectral assumption on the forcing modes, which ensures that the 
forcing is non-resonant, or in other words that the boundary layers remain 
stable. 

(ii) The above theorem holds for all values of the ratio v/e, but the as- 
ymptotic behaviour of ul depends on the scaling ofv/e. 

Note that, in the case when e>v, the effects of the boundary terms, even 
damped by the penalization, remain localized in the vicinity of the surface 
and thus do not contribute to the mean motion. 

If vje — * 00, the vertical dissipation damped by the penalization induces 
a strong relaxation mechanism, so that we expect the solution to be well 
approximated, outside from some initial layer, by a "stationary" solution to 
the wind-driven system. That initial layer should be of size O and the 

relaxation should be governed by the Ekman dissipation process (|2.5p . 



(2.6) 




in L-(R+,L 2 (T 2 x[0,l])). 
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(iii) If the forcing a bears on resonant modes only, then we are able to 
prove a global result. Precisely, assume that 

a{r, x h ) = <7+e ir (l, i) + a~e^ T (l, -i). 

Then there exists some destabilization profile v v solution of the heat equation 
113. 10\) such that 



(2.7) u e {t) 



<>xp [--L )K"(/) + <V(M 







m L°°(R+,L 2 (T 2 x [0, 1])) n L 2 (R+, L 2 (T 2 x [0,1])). 
In particular, for large times, 

u(t) « exp Mt) = 0[fl). 

Since (3 may be very large (see (|4.17p ), there is a destabilization of the whole 
fluid inside the domain as t — » oo. Note that the two convergences (12.6|) and 
(|2.7p are compatible, since with assumption (|4.17p . 



Vv = 0(v 3/4 (3) = o(l) in L 2 ([0, T] x T 2 x [0, 1]) 
for any finite time T > 0. 

2.4. Method of proof. 

Let us now give some details about our method of proof. As the evolution 
equation is linear, we will use some superposition principle, meaning that 
we will deal separately with the forcing and with the initial condition. 

• More precisely, we will consider on the one hand the wind-driven system 

d t u + ~P(e3 A u)AhU — vd zz u = 0, 
e 

V • u = 0, 

( 2 - 8 ) u lt=0 = 0, 

u\ z=0 = 0, 
U3|z=i = 0> dzUh\z=i = 0cr e . 

For that system, we will construct an approximate solution constituted of 
a boundary term u BL ' 1 localized near the surface, and some interior term 
v int,l^ which accounts for the fact that the vertical component of u BL ' 1 does 



not match the no-flux boundary condition at the surface (see Remark [2? 
("))• 

The convergence of the modes such that ^ 1 is then proved using a 
somewhat soft argument, which can be applied with a crude approximation. 

Concerning the quasi-resonant modes, for which = 1 and fc/, / 0, the 
situation is more complicated, and we have to build several correctors before 
reaching the adequate order of approximation. 
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• On the other hand, we will study the initial value problem 



d t u + -P(e 3 An) — ^A fe « - vd zz u 



V • u 



(2.9) 



u\t=o 

U\z=0 

u 3 \ z =i = 0, d z u h \ z=1 



1, 
0.. 
0. 



Here we will use, following [3], an energy method which requires to obtain 
a very precise approximation. A quantitative result about the required pre- 
cision is given in the stopping condition in the Appendix (Lemma [1]): when 
the approximate solution satisfies the hypotheses of LemmaHJ we put an end 
to the construction of correctors and conclude thanks to an energy estimate, 
whence the name 'stopping Lemma'. The approximate solution is actually 
obtained as the sum of two interior terms u mt that we seek in the form 



coming from the analysis of the linear penalization as an operator of L 2 , and 
two boundary terms u BL '°. We emphasize that in the case v = O(e), the 
construction of an approximate solution for system (|2.9p has already been 
dealt with by several authors (see [3J OS]); we recall it here for the reader's 
convenience, and further extend it to the case when v 3> e. 

Of course, in the nonlinear case the superposition principle does not hold 
anymore, and both systems (|2.8j) and (|2,9p will be coupled. 

The next sections are devoted to the proofs of Theorems 12.21 and 12.51 
We start with a precise description of the boundary layer operator B in 
Section We then build, in Section 0J the approximation and prove the 
convergence for the (possibly resonant) wind-driven system ()2.8|) . For the 
sake of completeness, we finally study the system (|2.9p which has already 
been dealt with in a number of mathematical papers. Let us recall that in 
both cases we need a refined approximation with many orders. We have 
then to iterate some process giving the successive correctors. Note however 
that we are not able to really obtain an asymptotic expansion leading to a 
more accurate approximation (in L 2 sense). At each step of the process the 
order of the resonances involved in the estimates is indeed increased, so that 
it is not possible to obtain convergent series. For more precisions regarding 
that point, we refer to the proof in Section 



u' 



,int 



°i N i e 



3. The boundary layer operator 
This section is devoted to the proof of Theorem 12.21 
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3.1. Non-resonant case. 

We recall that the boundary conditions are given by (|2.ip . and that we 
seek the boundary terms as a sum of oscillating modes, rapidly decaying in z. 
Our goal in this paragraph is to characterize these modes, or in other words 
to describe the propagation with respect to z of the boundary conditions 

W| z =o = $1, d z v h \ z=1 = 8\ . 

We will use the following Ansatz 

v(t,x) = v°(t,x) + v x {t,x) 

with 

v 3 {t,x) = ^2v J (n,k h ;x)exp ii-fi 

where \x and kh are the oscillation period and horizontal Fourier mode. 
We further seek V°(fi, kh) and V 1 ^, kh) in the form 

V°(n, k h ; x) = v°(n, k h ) exp(ik h ■ x h ) exp f -A(/x, k h ) 



(3-1) ' (J _ 

V 1 (n,kh,;x) =v l (ix,kh)exp{ikh -x h )exp ( -X(fj,,kh s 



feu 

so that they are expected to be localized in a neighbourhood of size 0{yfeD) 
respectively near the bottom and near the surface. Note in particular that, 
with such a choice, v° (resp. v 1 ) introduces only exponentially small error 
terms on the surface (resp. at the bottom). 

Plugging this Ansatz in the system (|1.4p we get actually 

\2~ . j,2 ~ * , hk 2 h - kfv 2 

ifivi — X vi + ek h vi -v 2 + eu — — = 0, 

A - euk h 

(3-2) . A . 2 . . ,2- , * , -fafah + k%vx 

ijj,v 2 - X v 2 + ek h v 2 + vi + eu -r^ — = 0, 

A - euk h 



\fev{ik\V\ + ik 2 v 2 ) ± Xv^ = . 

which expresses the balance between the forcing, the viscosity, the Coriolis 
force and the pressure. 

Denote by A\ the matrix corresponding to (13. 2D 



(in A 2 i rk 2 I euhk2 1 eVkl 

+ €kh + X 2 -evk\ ~ l ~ X 2 -evk\ 

eukn . , 2 ,2 euk\k 2 



i + TtT 'L.2 i^-X 2 +ek 2 h 



\ A 2 - evk\ h A 2 - euk 2 J 

Classical results on boundary layers are then based on the fact that \n\ ^ 1, 
which ensures that the matrix 

[i i 
-i fj, 
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is hyperbolic in the sense of dynamical systems, i.e. that its eigenvalues 
have non zero real parts. In particular, there exist two complex numbers 
A = A(/x, kh) with nonnegative real parts such that det A\ = 0. 

This feature, as well as general properties of the system, is therefore 
stable by small perturbation. The method consists then in neglecting the 
perturbation, i.e. the pressure and horizontal viscosity terms and to compute 
a solution to 

d t v + e 3 A v - vd zz v = 
with suitable boundary conditions. 



Now, if \/i | = 1, the matrix 



fx i 
-i ji 



admits zero as an eigenvalue, and we expect its behaviour to be very sensitive 
to perturbations. Actually we will distinguish between two cases 

• either kh ^ (0, 0) and we will prove that the same type of behaviour 
as previously occurs, with the difference that the decay rate A of the 
singular component is anomalously small. We will thus develop a 
general method, which can be used independently of the size of A 
(the classical method fails since the error depends on 1/A 2 ). 

• or kh, = (0, 0) and we have a bifurcation. The solution v is not 
localized anymore. 

Case when kh / (0,0). 

• Let us first introduce some notations in order to define an abstract 
framework to deal with. For the sake of simplicity, we omit here all the 
parameters /x and kh- 

Let A be such that det(A\) = 0, then there exists w\ such that 

(3.3) A x w x = 0. 

In other words the vector fields W® and Wl defined by 

/ _W\ \ t ( Z 

W x (t,x) = yfev ex.p(ik h ■ x h ) exp(i[i-) exp -A-= 

\ — 7-ik h - nix e V V ei/ 

(3.4) / A \ 

w \( t i x )=\ fl h WX I exp(ik h • x h ) exp(i/i-) exp ( -A ^ ^) ) 
\-^ikh-w x J e V J 

are exact solutions to (|1.4p satisfying respectively the horizontal boundary 
condition 

W x,h\z=o = W ^ exp(iA; ft , • x h ) exp(i/i-), 



A s ,. t s f\ 



d ^ W \,h\z=i = T^^A exp(ik h ■ x h ) exp(i(i-) exp 
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and 



d * w \,h\z=i = w x ex.p(ik h ■ x h ) exp(i/x-), 



W X,h\z=o = ^~ Wx ex P( ik h ' x h) exp(ifi^) exp (--^==\ 



We have moreover the following estimates (provided that -j= 3> 1 



(3.5) 

wl = o(i) L oo {R+ ^ m , w° x = o ((^Y /4 ) 

\\A / / L°°(R+,L 2 (n)) 



w^oim 1 ' 2 ) , wi=o m m 



1 L°°(R+,L°°(Q)) \\A J / L°°(R+,L 2 (n)) 



We intend to build one particular solution to (|1.4p satisfying the horizontal 
boundary condition 

v h \ z=0 = o h , 
d z v h \ z=1 = 8\. 

Hence, we only have to find (for all /j, and kh) some w x - and w x + constituting 
a basis of C 2 . 

• In order to determine some suitable w x - and w x + , we have to get some 
asymptotic expansions of the eigenvalues and eigenvectors of A x ({i,kh). 
In view of the previous paragraph, at leading order, we have 

so that 

det(A A ) = (i/z - A 2 ) 2 + 1 + o(l) = 
for (A~) 2 = i(n + 1) + o(l) or (A+) 2 = i(fi - 1) + o(l). We further have 
u> A - = (1, — i) + o(l) and u> A + = (1, i) + o(l) 

For ^ 1, we choose A - and A + to be the roots of det(^4 A ) = with non- 

: previous asymptotic equi 

det(iu A - , w x +) = 2i + o(l) 



negative real parts. The previous asymptotic equivalences are then enough 
to prove that 



from which we deduce that (w x - ,w x +) is a (quasi-orthogonal) basis of C 2 , 
and that we have uniform bounds (with respect to e sufficiently small and 
v bounded) on the transition matrix P and its inverse. 
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For fj, = 1 we expect A - to be given by (A~) 2 = 2i + vj~ with rj_ = o(l), 
and A + to be given by (A + ) 2 = rj + with rj + = o(l) 

_ («„ - A 2 + £ * 2 + ^) (<, - A 2 + ^ - ^ ) 



}? — evk\j \ }? — tvk\ 
= [-i-V~+ ekl + ^1 f-i - rf + dfcg - ^ 



2i J \ ' 2i 

evh 

2i J V ' ~2i 



, 2 euk 1 k 2 \ ( . + 2 e^i&2 



and 

^ ( 1 + ^ l ) +0 ( eV ^ + ) 2 ) 

from which we deduce that 

rf = ek\ + -eukl + o(e) 
= efch + |^ + o(v^) + o(e). 

We have then 

(A~) 2 = 2i + 0(e). 

On the other hand, a discussion taking into account the relative sizes of e 
and v shows that 

(A+) 2 ~ ek 2 h if v « e, (A+) 2 ~ ±^\k h \(l + *) 
while an easy argument of homogeneity gives 

(A+) 2 ~C(^)eif^~e, 

for some constant C(kh), depending only on k^. Thus there exists a constant 
C(kh) such that 

c(k h y l < \\-{i,k h )\ <c(k h ), 

C(k h )-\e + V^) 1 ' 2 < |A + (1, k h )\ < C(k h )(e + y^) 1/2 . 
Plugging these expansions in the formula of A\ leads then to 

w x - = (l,-t + 0(e)), 
w x+ = (l,i + 0(^) + 0(e)) 

In particular we have 

det(w A -,u; A +) = 2i + O(e) + 0(y/eu) 
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from which we deduce that (w\- ,w\+) is a (quasi-orthogonal) basis of C 2 , 
and that we have uniform bounds (with respect to e and v sufficiently small) 
on the transition matrix P and its inverse. 

For fi = — 1 we have in the same way 

C(k h )-He + V^) l/2 < |A-(-l, k h )\ < C(k h )(e + 

C(k h y l <\\+(-l,k h )\<C(k h ) 

and 

w x - = (1, -i + O(v^) + O(e)), w x+ = + 0(e)) 

from which we deduce uniform bounds (with respect to e and v sufficiently 
small) on the transition matrix P and its inverse. 

• We then define V°(/i, kh) and V 1 ^, kh) by 

(3.8) V j {fi,k h ;x)exp (ifJ^j = oP_W 3 x _(t, x) + a{W J x+ (t, x) 

where W 3 X is defined in terms of w\ by (|3.4p and the coefficients oP_ and oP + 
are defined by 

(3.9) (d,a> + ) = p- 1 Si(n,k h ). 

Case when fc^ = (0,0). 

That case is strongly different since there is no term of higher order in 
432} : 

ifx-X 2 -1 
1 ifi- X 2 

For ^ 1 we use exactly the same arguments as previously and define 
v j (^,0) by formulas (ETg|l (|3T9l) . 

Remark 3.1. When \fj,\ = 1 we cannot find a basis of eigenvectors (w\- , w x +) 
with 5R(A _ ) > and 5R(A + ) > 0. One of the eigenvalue is necessarily 0, and 
thus the corresponding solution has no decay in z. In other words we do 
not expect the boundary terms to be localized in the vicinity of the boundary 
uniformly in time. 

The assumption that the boundary condition is non resonant ensures 
however that there is no such contribution. 
If |//| = 1 we have 

= 2(j,i and = 

with 

w x - = (1, — *) and w x + = (1, i) . 
If we define as previously by (|3.4|) . and a± by (|3.9p . we have 



16 ANNE-LAURE DALIBARD AND LAURE SAINT-RAYMOND 

Setting, for j = or j = 1, 

V j (/J,,0;x)exp = ^-^{-^ (*> x ) 

we can check that it is an exact solution to (|1.4I) . which further satisfies the 
required horizontal boundary condition. 

3.2. Resonant case. 

Let us then focus on the resonant part of the motion. The singular com- 
ponent u etres of the velocity is a 2D vector field (depending only on t and 
z), so that (jl.4p can be rewritten 

d t u e 

,res ,res 

A e 3 - vd zz u e res = 0, 

e 

meaning that the pressure is constant. 

• Therefore the equation can be filtered by a simple change of unknown : 






A straightforward computation leads then to 
(3.10) d t v v - vd zz v v = 0, 

which is nothing else than the heat equation with small conductivity v. We 
therefore expect the boundary effects to remain localized (in L 2 sense) in 
layers of size 0{yfu£) near the boundaries. 

• Let us then introduce a boundary layer approximation 

VL,res = v L,res + v L,res 

for v u . The heat equation on v 3 Lres is supplemented with the boundary 
condition 

ii° — A ft ii 1 — A 1 

u L,res\z=0 ~ u L,resi u z u L,res\z=l ~ u L,res 

and the initial condition 

Vl U n = . 

L,res\t=0 

Notice that once again we do not enforce boundary conditions on both sides 

J 

L,res 

L,res 



for v J L res : the trace of v 3 at z = 1 — j will be imposed by the exponential 
profile condition. We indeed seek vi in the form of self similar profiles 



(3.11) V° L tres = if° [-=), d z v\res = V l 



z \ „ i 1 /(l - z 



vt J V yvt 

We then get 

-\xjix) - V "{X) = 0, 
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from which we deduce that 



^(X) = ^(0)exp(-ix 2 ) , 

-^ + °V(0)exp (~yA dY. 



and 

<P(X) = 
We thus choose 

Note that, in order that v\ res satisfies the heat equation (|3.10p . we have to 
further impose that v\ Tes {— do) = 0. 

• We deduce that 

1 f + °° Y 2 

V L,res = -/^ S L,res I &"~dY 



~** (^) exp B (^) 



Similarly, we have 



with 

991 (w) ^ u (w) exp H 

Therefore 

1 / x 4(^) 3 /2 / 1/1 

Z ) ^.resC 1 - ^) exp I -- I - 



. 2^ 

Z 



.'vt 

In particular v J L res is exponentially small outside from a layer of size 0(y/ui). 

3.3. Continuity estimates. 

We now turn to the derivation of the estimates of Theorem 12.21 
Thanks to the previous paragraph, the resonant part of the boundary 
layer, namely v 3 res defined by (|3.11|) . satisfies the third estimate in f|2.3[) . 
We then split v 3 — v 3 res according to the size of the boundary layers 
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By definition of ai(fi, kh) and W{, 



, we then obtain the estimates 



-J 




for the classical boundary layer, and 




y-e + jeu) 1 / 4 



e + y/eu 



) 



4 



\K 



for the quasi-resonant boundary layer. 



4. Study of the wind-driven part of the motion 



This section is devoted to the proof of Theorem 12.51 in the case where 
the initial data 7 vanishes. In other words, we study here the asymptotic 
behaviour of the system (|2.8p . Our goal is to prove that under a technical 
scaling assumption which will be precised later on, the solution u of (|2.8p 
converges towards zero in L^ C (R + , L 2 {uj)) as e, v — > 0. 

As explained in section [21 the method of proof relies on the construction 
of an approximate solution u app , defined as the sum of boundary layer terms 
obtained thanks to Theorem l2.2l and interior terms which will be determined 
by a filtering process. The presence of these interior terms is due to the 
fact that the vertical components of the boundary layer terms constructed 
in Theorem 12.21 do not vanish on z = and z = 1. More importantly, 
the traces of these boundary layer terms do not satisfy the assumptions 
of the stopping Lemma [U in Appendix B, which quantifies the order of 
approximation required for u app . Hence in general, the approximate solution 
is constituted of several correctors, which all vanish in I? norm. 

The different modes of the wind stress a will be treated independently of 
each other. Indeed, in the case where the stress a does not have any quasi- 
resonant mode, it will be sufficient to construct a very crude approximation, 
constituted merely of one boundary layer term and one additional corrector. 
On the other hand, the vertical components of the quasi-resonant boundary 
layer terms have a much larger trace on z = 1 and z = than the classical 
ones, as can be seen in inequalities (|2.3p . Consequently, the quasi-resonant 
part of the stress a will require a much more refined approximation, with 
several orders of boundary layer terms and interior terms. 

The organization of this section is as follows: first, we give in paragraph 
I4.1l a general convergence result for the system (12.8|) . Then, in paragraph l4.21 
we construct the first orders of the approximate solution u app . In paragraph 
14.31 we conclude in the case when there is no quasi-resonant mode \fi\ = 1, 
kh 7^ 0. At last, we prove the theorem for the quasi-resonant part of the 
stress a in paragraph 14.41 At each step, we give some sufficient assumptions 
on the parameter 0, and at the end of the proof, we only keep the most 
restrictive ones, which will lead to the scaling assumption (14. lTj) . 
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4.1. Some stability inequality for the wind-driven system ( 12.81) . As 

mentioned in Section 2, for the non-resonant part of wind-driven system 
(|2.8j) . we will only need a rather crude approximation of the solution. We 
have indeed the following 

Proposition 4.1. Denote by u e the solution to \2. 8\) and by u app any ap- 
proximate solution in the sense that 

dtUapp + ^P(e3 A U app ) - A h U app - vd zz Uapp = TJ, 



(4.1) 



^ ' Uapp — 0, 
u app\t=0 = Vinii 

u app,3\z=o = 0, u a pp,h\z=o = erjo, 

U app ,3\z=l = 0) dzU app ,h\z=l = + Vl, 



with r] — ► in L 2 ([0, T] x u), rji n i — > in L 2 (uj) and rjo, u 3 ^r]i, edtrjo —* in 
L 2 ([0,T] x uj h ). Thenase,u^0, 

IK - u app\\L°°([0,T],L 2 (uj)) —> 0, 
||Vfc(u £ - U a p p )\\ L 2( [0yT]xLu) + y/v\\d x (u e - U a pp)\\ L 2([0,T]xuj) ~> 0. 

• The first 
Uapp + w satisfies 



Proof. • The first step consists in building a family w such that v app 



dtVapp + ^IP(e3 A v app ) - A h v app - vd zz v aV p = C> 

^ ' V a pp — 0, 
v app|i=0 = Cinii 
v app,3\z=0 = 0] v app,h|z=l = 0) 

v a pp,3\z=i = 0, d z v apP: h\z=i = (3o- e + 771. 

with Qni -> in L 2 (a;) and C ^ in L 2 ([0,T] x oj h ). 

In order to do so, we just apply Lemma Q] in the Appendix with 

£0 = -er/o, 6% = and = 0. 

A simple computation allows then to establish all the properties (|4.2p . 

• The convergence is then obtained by a standard energy estimate. Com- 
bining (|4,2p and (|2.8p . and integrating by parts lead indeed to 

-||(«e - Vapp){t)\\ 2 L 2 + / [|Vfc(K - V a pp)\\ 2 L ads + U [ \\ Vh{v>e ~ V app ) \\ 2 L ads 

1 Jo Jo 

1 2 t 



+^ / l|(«e - V ap p)h\ z=l(t) \\L 2 (ui h ) IK \\L 2 (uj h )ds 
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• To conclude we therefore need to estimate the trace (u t — v a pp)h\z=i in 
L 2 {ijjh) in terms of the H 1 norm of u e — v app . By Sobolev embeddings and 
the Cauchy-Schwarz inequality, we have 

v 1,2 \\{u e - v app ) h \ z=l f L 2 {uih) < C\\u e - v a p P \\ 2 L 2 {uj) + v\\d z (u e - v ap p)\\ 2 L 2 H 

Plugging that estimate in the previous energy inequality, we get 

2 II («e -«ORp)(t) Ilia + ^ II Vh(n e -V app ) ||i 2 ds+ - jf ||«9 2 (« e - U app ) ||^ 2 ds 



<^IKmi|li 2H + i I \\<\\h H ds + v 3/2 I \\m\\h M ds 



1 * 



JO 

using again the Cauchy-Schwarz inequality. We conclude by Gronwall's 
lemma 

r-t rt 



app)(t)\\ 2 L 2 + / \\V h (u e - Vapp) Wtfds + V I \\d z (ll € - V a pp)\\ 2 L 2ds 

JO Jo 



2"l |1 '"' ~ '"' 

e 2Ct „ „, 1 ''' 

'o jo 



< VlWI^) + 5 /„ llClli 2 He 2C(i - s) ^ + ^ 2 jf IMi^e 20 ^* 



which proves that M e — t> ap p converges to in L^ c (R + ,L 2 (a;)). Theorem 



will be proved in the case when 7 = if we are able to build some 
approximate solution u app that converges strongly to as e, v — > 0. □ 

Remark 4.2. The above proposition can be slightly modified if one wishes 
to work with a source term n belonging to L 2 ([0,T], H~ 1 (lo)), for instance. 
In this case, following exactly the same argument as in the proof above, the 
relevant assumption on r/ is 

( 4 -3) -?=\\v\\l^([o,t],h-^(u>)) = ase,v-*0. 

4.2. The first order terms of the approximate solution. In order to 
obtain some approximate solution to (I2.8p (in the sense (|4.ip of the pre- 
vious paragraph), we will essentially need to construct the boundary layer 
term and some small corrector to account for the vertical component of the 
boundary condition. 



• We define with the notations of Proposition | 

U ■ =D(U,pa)=U ' + U ' + u res ' . 

Since we assume that a has a finite number of horizontal Fourier modes kh 
and of oscillating modes /i, by Lemma 12.21 we have 
11 -Bi.l 1 



li» < C\\a\\ L 2 {u)h) f3{evf/\ 
~hi L,1 \\i?{w) < C\\a\\ L 2^/3(eu) 5/i , 
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and for the quasi-resonant modes with |//| = 1, kh ^ 0, 

/ \ 3/4 

Wu^Wl^) < C\\a\\ LHwh) f3 < C^/\ 

5/4 



As for the resonant modes = 1, kh = 0, we get 

''/iIIl 2 ([0,T]xu,) < C/3T' 5/4 ||cr|| X ao([ 0)T ] !l/ 2( a , ;i) )i/ 3/4 , 



Hence u BL > 1 vanishes provided 
(4.4) 0u 3/4 = o(l) ase,i/->0. 

Furthermore, using the explicit formula for B, we get 

ll^ijifsK) = 0(/3(ev)) and II^w^Jh^^) = 0(/3i/) 
HfiJJJjII-ff'K) = 0(p(ev) N ) and ||ftt^f'i ||ir»( Wh ) = 0(/3(ev) N ), 

ll^ij^K) = 0(/?(^) 1/2 ) and ll^u^ij^^) = O (z 9 ^) 
HwJSillff'K) = 0{p(ev) N ) and ll^u^JI^^) = 0{p[ev) N ) 



for any integer iV, and uniformly in time. As a consequence, u^^, u 3 i ^ Q 
and u^,^ satisfy the conditions of the stopping Lemma [T] in the Appendix 
as soon as /3v = o(l), which is always ensured by hypothesis (|4.4|) . We 
denote by w the function defined in Lemma [T] with 



5° h = 0, 5 h = 

A _ -BL,1 ~BL,1 A l _ -BL,1 
°s - u 3\z=0 U 3\z=0> °3 - u 3| 2 =r 



• The term uf, L '\ , on the other hand, does not match the conditions of 

3|z=l> ' 

Lemma [TJ We therefore introduce some corrector v mt ^ to restore the zero- 
flux condition. We first define its vertical component 

int,l _ _ ~BL,1 
V 3 - U 3\z=l Z ' 

then its horizontal component in order that the divergence-free condition is 
satisfied 



V7 f A \ — l~BL,X 

V h (A h ) %|z i r 
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Note that for kh = 0, v 1 '^' 1 is identically zero. In any case, we get easily 
that 

IK M IIl°°([o,oc),^m) = O h3 e+ e ^_ ) = 0(/3v), 



\\d t v in > ||l°°([o,oo),l 2 (^)) — O P\ - 



With the above notations, the first order of the approximate solution is given 
by 



u l pp = u BL,l + w + v int,l_ 

4.3. Proof of Theorem 12.51 when there is no quasi-resonant mode. 

If there is no quasi-resonant mode (see the precise definition in the previous 
section), namely if u BL ' 1 = 0, we then claim that u\ pp satisfies the required 
conditions. We indeed have clearly 

U app,3\z=0 = U app,3\z=l = ® 

by definition of w. Notice that in this case v m1,1 = 0. We further have 

d zu\ ppMz=1 - fo e = 0, 

ll«Jpp,fc|*=olU a («jO = 0(l3(eu) N ) and lldf^^olll^) = 0(/3(eu) N ) 

for all N . We also have for all t > 

hlpp(t)\\L^) < \\u BL ^{t)\\ L , {u}) + \\w(t)\\ L2{w) = 0(f3{euf/ A ) = o(l). 

It remains then to check that the evolution equation is approximately 
satisfied. We have 

dtulpp + ^P(e 3 A u\ pp ) - A h ul pp - vd zz u\ pp = 0{v(3) L 2 {[ ^ T]xu}) = o(l) 
supplemented with some initial condition 

u 1 applt=0 = O((3(ev) 3 ^) L 2 H . 

We therefore apply Proposition 14.11 and conclude that u^ pp has the same 
asymptotic behaviour as the solution of 

d t u + -P(e3 An) — AhU — vd zz u = 0, 

V-« = 0, 

n 3 | 2=0 = 0, u h \ z=Q = 0, 
n 3 | z=1 = 0, d z u h \ z=1 = /3a € . 

Since u^ pp vanishes in L^ C (R + , L 2 (lu)), Theorem 12.51 is proved when 7 = 
and when there is no quasi-resonant mode in the forcing a. 
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4.4. Proof of the Theorem in the quasi-resonant case. For the quasi- 
resonant modes \fj,\ = 1, the influence of the forcing is much more extended 
inside the domain. In particular, the defect 

S =d t v int ' 1 + -e 3 A v int ' 1 - AhV™*' 1 
e 



M=±l k h 



-v 2 {n,k h ,z) 



+ 7 EE vt'^M 

M=±l k h \ / 

does not converge strongly to in I? norm. It is however expected to have 
rapid oscillations, and thus to converge weakly to 0. The standard method 
to deal with such a problem consists then in building some corrector which 
will be small in 1? norm in contrast with its time derivative which has to 
compensate the previous defect. 

More precisely we will use the small divisor estimate stated in Appendix B. 
For K > arbitrary, denote by Su^ 1,1 = Y2i wie~ l e Xl Ni the solution to 

fySu**' 1 + -F(e 3 A 5uf' 1 ) - A h 6^ - vd z M^ = -P*(£), 
supplemented with the initial condition 

Su T\to = °- 

The notation stands for the projection onto the vector space generated by 
{Ni, \ l\ < K}. The idea is the to choose carefully the truncation parameter 
K, depending on e and u, so that both Su^' and the error term P(S) — 
Pi^(S) are small in suitable Sobolev norms as e and v vanish. 

• Let us first derive the equation on w\. For |/| < K, wi is the solution of 

dm + llrfm + v'lhfwi = -e iA ^(iV z |£) 
where v' = ir 2 u. Direct computations give for 1^ ^ 0, fj, = ±1, 



v^ivJ^z) = 5 3 (n,l h )z, 

where 

=il3{ev) —- , 

where a* and w\ were defined in the previous section by (|3.9j) and (|3.3p 
respectively. Notice moreover that A M satisfies the estimates (|3.6j) - fj3.7|) . so 
that in general, 

(A^)(^,A;,)- 2 = 0(M- 1 / 2 ). 
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Moreover, 



(4.6) 



N, 





*-fe£-(-i) /s wo 



il h -x h 



2vr 



if h + 0, 
else. 



We thus have 



(4.7) = I £ 



2tt 



d t wi + {\l h ? + v'\h\ 2 )m 



ttZZ, 



• We now estimate the different terms and explain how to choose the trun- 
cation parameter K. Notice first that by truncating the large frequencies in 
Z, we have introduced a source term in the equation. Precisely, du 1 ^' 1 -\-v mt ' 1 
is a solution of equation (jl.4p with a source term equal to 



(S - PS) + (PS 



K 



S). 



The term S — PS belongs to by definition of P, and thus for all u G Vo, 
we have 

I (S — PS) ■ u = 0. 

J u 

As for the remainder term PS — P^S, we have 

||PS - PxS|| i o O (( 0iOO ) iL 2 (a ,)) < C(3^K~ 3/2 , 

||PS - PifS|| L oo(( 0)Oo)iH -i( w) ) < Cf5^~-K~^l 2 . 

With a view to apply Proposition 14.11 or its variant sketched in Remark 
4.21 we need the source term PS — P/fO" to be either o(l) in 1? norm or 
o(yfv) in H~ l norm as e, v — > (see condition (|4.3p ). Precisely, according 
to Proposition 14.11 and Remark 14.21 the parameter K should satisfy either 



(4.8) 



/3^V 3/2 = o(l) 



as e, v —* , 



(4.9) 
or 

(4.10) -L 0* f^K-V 2 = = o(1) ag y ^ 

V e v e 

On the other hand, we apply Lemma [2] to get 
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For further purposes, we have to choose K such that the H s norm of du 1 ^' 1 
satisfies 

J^W^K^WH'iw) = as e,i/ ->■(), 

and such that at least either (14.91) or (|4.10j) is satisfied. We distinguish 
between the cases when v is large (say v > e) and z/ is small (say 1/ < e), 
which yield different values for K. 
- If v < e, we choose X so that 

for some s > 3/2, which yields 



K = (eu) 2 ( s+2 > . 

With this choice, we have 

Now, assume that satisfies the following assumption 

3(a ,ai) £ (0,oo) 2 , a < 5/7 and ax > 2/7, 3C > 0, 

< e =>■ /3 < Ci/ _ao e ai . 

We choose sq > 3/2 such that 

s o + \ 
2(s + 2) 

s o + I 

"1 ; \" > 0, 

2( So + 2) 

and we have, as e, v — > 0, 

(4.12) ||PS - P^S|| i2 + 7?||<5«S t,1 ||^ (a ,)=o(l). 

- Else, we choose K so that 

for some s > 3/2, which yields 

i_ 

K = [Vyfe) 

Assume now that (3 satisfies the following assumption 

3(q!q, a%) e (0,oo) 2 , a < 5/9 and ax > 2/9, 3C > 0, 
^ 4 ' 13 ^ v > e (3 < Cv~ a °e a \ 
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We then choose so > 3/2 such that 

s o + \ 
2(s + 3) 

1 - ao - | > 0, 
s + 3 

and we have, as e, f — > 0, 
(4.14) 

— ||PS - PxS||loo((o,oo),H- 1 (o;)) + ■y^~ll^ n A' < ' 1 |l^ oo ((0,oo),H s 0(w)) = f(l). 

We emphasize that this method remains valid when v = O(e); however, 
if v = e, condition (|4.13j) is more restrictive than (14.111) . 

• Because of the terms v int ' 1 and Su^' 1 , the horizontal boundary condi- 
tions are no longer satisfied at z = (notice however that they are satisfied 
at z = 1). Thus, we construct another boundary layer term, which we 
denote by 5u BL,1 : such that 

c BL,1 jo/ int,l c mt.l n \ 

5u =S(-^ |z L -<5n K ^ =0 ,0). 
The above definition is not entirely licit, since _ takes the form 

[i|<ir v 7 

where the vector is defined in Appendix A (see (|6.7p - (|6.8p . Hence 

^KA|2=o depends on the fast time variable t/e, but also on the slow time 
variable t through the coefficient ihi. In the definition of 5u BL,l : we forget 
the time dependance of wi, and consider the coefficients rbi as constants. 
Consequently, the boundary layer term Su BL}1 is not an exact solution of 
equation (jl.4p . but there is an error term depending on dtwi. This error 
term will be estimated later on. 

We now decompose <5it' Bi ' 1 into 5u BL,x = 5u BL)l + du 31 "' 1 as in Lemma 
12. 21 the term 5u BL)1 is due to the modes kh ^ 0, = 1, and thus depends 
only on v mt ' 1 , since |A&| < 1 if kh ^ 0. Notice that there is no term Ju^s' 1 
because v mt > l (n, lh, z) = for fc^ = 0, wi = for kh = 0. 

According to the estimates (|2,3p . and provided (|4.13p holds, we have, for 
all t > 0, 

[|*£ L4 (9[|l»( u ) < CII^WII^^) (ci/)*+ qifti^COII^Cci/) 1 ^ 

< CV) 1/8 , 

ll*uf x,1 (*)llx» < IK^WIbn M 1 +c||^ 1 (t)||^ 0+ iM 3/4 

< C(ei/) 3/8 + C(e^) 3/4 (i/v^)"^ 
Thus 5u BL > 1 vanishes in L°°([0, T], L 2 (uj)). 



j 

1+3 
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Let us now estimate the error term in equation (11. 4p due to the time 
dependance of wi. According to (|4.7p . there exists a constant C such that 



\9 t Wl\ < 



so that 5u BL ' 1 is an approximate solution of equation (|1.4[) . with an error 
term which is bounded from above in L 2 ([0, T] x co) by 

Hence, the new condition on [3 is 

(4.15) = o (V 3 / 4 e 1/4 ) as e, v -> 0. 

Notice that (|4.15p immediately entails (|4.4p . 

• Let us now check that the remaining boundary terms are all sufficiently 
small to conclude. To begin with, the terms 5u B h^ v 8u^=i are exponen- 
tially small, and thus satisfy the hypotheses of Proposition 14.11 and Lemma 
[1] respectively. We now prove that under conditions (|4,lip - (|4.13p . 5u^= 
also satisfies the assumptions of Lemma [TJ Using the construction of the 
previous section, it can be checked that 5u BL ' 1 is given by 



5u BL '\t) 



E E E e-^al(-X k ,k h )W^ 
\k h \<N \k 3 \<K 







+ E E j^-MkhWz 

\k h \<N ^{-1,1} 
where the coefficients satisfy 

VfcGZ 3 , \a^(-X k ,k h )\<C\w k \, 
and \a°_^(n,k h )\ <C 5 3 (n,k h ) . 

Recalling the expression of W® (see (|3.4p ). we infer that for all t, 



< 



c E E E i^Itw-x, 

+ C^) Yl E 

\k h \<N 

< CV^ Yl E \k\\w k \ + Cf3ev, 

\k h \<N\k 3 \<K 
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and thus, using the Cauchy-Schwarz inequality (recall that so > 3/2 and 
that N is bounded) 

1/2 



5u 



BL,l 
3|z=0 



< CJeu sup 
L«([0,n,i a (wfc)) te[0,T] 



.lfc 3 l<tf 



+ C(3ev 



r intA 

8u K 



+ Cfiev. 



L°°([0,T],H s o) 

Hence, under conditions (|4.1ip - (|4.13p and by definition of K, the remaining 
boundary term 5u.^ z ^ Q satisfies the conditions of Lemma [TJ 



We now consider 5u 



~BL,1 
'3\z=Q- 



which is due to the modes fi = ±1, ^ in 



\5u 



BLA 



3| 2 =0lli°°([0,T],L 2 (a;O) 



< c 



eu 



(ei/)V4 + ^ 
< C/3(ei/) 3 / 4 . 



i n ,* n *>i ii 

\ v IIl°°([0,T],L 2 (o;)) 



3/2 



|o , ||L«([0,T],La(« fc )) 



Hence Su B ,^ satisfies the assumptions of Lemma [U provided (|4.15p is sat- 
isfied. 

Thus, we slightly modify the definition of the function w given by Lemma 
[U so that the boundary conditions are now 

s° h = o, 4 = o 



°3 



BLA 



A BLA , 
~ U 3\z=0 ~ 0n 3| 2 =0' °3 



-U 



BL,l 



Su 



BLA 



"3| 2 =1 UUj 3\z=l- 

• We then claim that under hypotheses (|4.11|) . (|4.13|) and (|4.15|) . 

u apP = u BL ' 1 +w + v int ' 1 + 5uf> 1 + 5u BL > 1 
satisfies the assumptions of Proposition 14.11 We indeed have clearly 

u app,3\z=0 = u app,3\z=l = 

by definition of v int ' and w. We further have, for all N > 

\\d z u amh \ z=1 - (3a e \\ L 2(u h ) = 0((eu) N ), 
\\uapp,h\z=o\\L2(u h ) = 0((eu) N ) and \\d t u apPth \ z=0 \\ L 2 {uJh) = 0((eu) N ). 
We also have for all t > 

\\u app (t)\\L*M < C^l i + 0{eu) 1 '\^eY^)=o{l). 

By definition of the different terms, the evolution equation is approxi- 
mately satisfied, up to an error term of order o(y / F) in L°°((0, oo), H 1 {uj)), 
and another one of order o(l) in L 2 ((0,T) x to). 
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We therefore apply the variant of Proposition 14.11 sketched in Remark 14.21 
and conclude that u app has the same asymptotic behaviour as the solution 
of 

d t u + ^P(e 3 An) - A h u - vd zz u = 0, 

V-tt = 0, 

( 4 - 16 ) u lt=0 = 0, 

u 3 \ z=0 = 0, u h \ z=Q = 0, 
«3|*=l = 0, d z ii h \ z=1 = f3a e . 

Thus the solution of (liTToT) vanishes in L°°([0, T],L 2 (o>)) norm as e, v — ► 
with (e,u,/3) satisfying KTTh . (14TT3J) and (l4T5j) . 

• We conclude this paragraph by giving a scaling assumption on (3 which 
entails all three conditions (14. lip . (|4.13p and (14.15p . Assume that the pa- 
rameter (5 is such that 

(4.17) 3a G (o,^j , f3 = O^e 1 / 4 ) ase,^ 0; 

we now check that each of the assumptions (14. lip . (I4.13P and (|4.15p are 
satisfied. 

First, it is obvious that 

I/ 3/4 e -l/4 j9 = O ( l/ 3/4-a ) = 

since 3/4 — ao>l/6>0. Hence (14.15P is satisfied. 

We now tackle condition (14, lift ; since ao < 7/12, there exists positive 
numbers (q!q,Q!^) such that 

a' < 5/7, a\ > 2/7, and a' — a\ = ao — — • 

In view of (|4.17|) . there exists a constant C such that 

j3 < Ci^" a o +a 'i-2e3 

< C v-< (^f^' 1 e a K 

Notice that a\ > 1/4, and thus if v < e, we deduce that 

Hence we have proved that <|4.17p =>- (|4.1ip . 

The treatment of (|4.13p is similar. We first choose positive numbers ot'^a" 
such that 

o.q < 5/9, 2/9 < a" < 1/4, and a' Q ' — a" = ocq — -. 
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Then if u > e, we have 



(3 < CV- Q o'+ Q 'M e 3 
< C i>-'<> [~)*~ a{ e"< 



Hence we also have (|4.17|) =^ (|4.13p . and eventually, we deduce that under hy- 
pothesis (|4.17|) . the solution of (|2.8|) converges towards zero in L°°([0, T],L 2 ) 
for all T > 0. 

5. Study of the dissipating part of the motion 

This section is dedicated to the rest of the proof of Theorem 12.51 Ac- 
cording to the preceding section, there remains to define the term u Dirichlet ; 
which is an approximate solution of (jl.4p . supplemented with the following 
boundary conditions 

Dirichlet n Dirichlet n 

U h\z=0 — U ' U 3\z=0 — U 

o Dirichlet n „, Dirichlet n 

°z U h\z=l - U ' M 3|z=l - U > 

„, Dirichlet 

u \t=o = 7- 

This point has already been investigated by several authors, see for in- 
stance [1]: the idea is to construct an interior term, denoted by u mt , which 
satisfies the evolution equation up to error terms which are o(l), and a 
boundary layer term, denoted by u BL , which restores the horizontal bound- 
ary conditions violated by the interior term. We emphasize that in order 
that the equation and the boundary conditions are satisfied up to sufficiently 
small error terms, we need to build some second order terms in both u mt 
and u BL . 

The organization of the section is as follows: in the spirit of Theorem 
12.21 and Definition 12.31 we first define an operator U, which allows us to 
construct an interior term, given arbitrary vertical boundary conditions. 
Then we explain how to choose the boundary conditions for the boundary 
layer term and the interior term in order to retrieve (jl.2p and (jl.3p with 
a = 0. In the last paragraph, we build one additional boundary layer term, 
and we prove Theorem 12.51 thanks to an energy estimate. 

Throughout this section, we use repeatedly the following norm: if 5 £ 
L°°([0,oo) x [0,oo),L 2 (u h )) is such that 

5(t,T,x h )= £a(-A fc A;t)e <fcfc - Xh e-* A * T , 
\k h \<Nk 3 ez 

where r stands for the fast time variable t/e, then 

ll*(V)H.:= I E E l^l 25 5(-A fc A;t)'° l 
\\k h \<Nk 3 ez 



ROTATING FLUIDS WITH RESONANT SURFACE STRESS 



31 



5.1. Construction of the operator U. Let 5\ and 5® in L°°([0, oo) x 
[0, oo), L 2 (u>h)) be such that 

(5.1) 4(t,r,x h ) = y, Y,%(- x ^ k ^ eikh ' Xhe ~ iXkT > 

\k h \<Nk 3 eZ 

and let 7 S Vo. In practice, the functions 5% and 5® will not be arbitrary, and 
will be dictated by the expression of the boundary layer operator constructed 
in the third section. In fact, we will see that 6% = 0, so that the expression 
of u int below is simpler, but we have preferred to keep an arbitrary value 
for 5g in order not to anticipate on this result. 
We define the operator U by 

where u mt is an approximate solution of equation (jl.4p and satisfies the 
following boundary conditions 

(5.2) u|f =1 = V^ 3 \ 

(5.3) n^ =0 = v^3°, 

(5.4) ufl = 7 + o(l). 

We emphasize that conditions (|5.2p - (|5.3|) will be satisfied exactly (without 
any error term). Of course the above conditions are not sufficient to define 
the term u mt unequivocally. We merely define here a particular solution of 
this system, which is sufficient for our purposes. 

The explicit construction of u int requires three steps: first, we exhibit a 
divergence-free vector field v mt, ° which satisfies the vertical boundary con- 
ditions (|5.2p - (|5.3|) . but not equation (|1.4|) . and then we define a function 
8u int, °, which satisfies homogeneous boundary conditions, and such that 

(5.5) u int := exp (-^\ uf l + 5u int '° + v int '° 

is an approximate solution of (jl.4p . supplemented with the initial condition 
(|5.4p . As usual in this type of problem, we first assume that exp {—t/e) u™ 1 
is the preponderant term in u int , and thus we begin by deriving an equation 
for the corrector term 5u int, ° involving u™'. Ultimately, this will allow us 
to write an equation for it™ 1 ". In the third step, we prove that the function 
§ u int,0 thus defined is of order 0{y/ve) in I? . 
• A natural choice for v mt '° is 

( ] I < nt '° = v^V.A, 1 [5° - 61} . 

(Note that v tnt, ° is not uniquely determined by (I5.2p - (|5.3I) ). We denote by 



V 



(/x, kh,t, z) the Fourier coefficient of v int, °, that is 

j int '°(t, x) = Y v int '°(v, k h ,t, z) exp(ik h ■ x h ) exp f i-ji ) 
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The fact that v™ ' u + means that a small amount of fluid, of order ^/eu5 J 3 , 
enters the domain (or the boundary layer, depending on the sign of the 
coefficient). This phenomenon is called Ekman suction and v™*' is called 
Ekman transpiration velocity. This velocity will be responsible for global 
circulation in the whole domain, of order (ev)^ , but not limited to the 
boundary layer. 

Furthermore the Ekman suction at the bottom has a very important effect 
in the energy balance. The order of magnitude of v J \ Vu BL \ 2 in the Ekman 
layer is indeed 0(-v/t) 3 so that the Ekman layer damps the interior motion, 
like a friction term. This phenomenon is called Ekman pumping. We 
therefore expect that the weak limit flow of (|1.4p in the high rotation limit 
is not determined by the formal equations (|1.6p but by a dissipative versions 
of this equation. 

• As in the previous section, we seek 



(5.7) 
(5.8) 
so that 



cxp 



-L ) u %* 



5u 



int.O 



£ S Cl (t)e- iX ^N l: 

26Z3 



d t + -L- vd z 
e 



A, 



exp 



L ) u\ nt + 5u int '° 



zez 3 



«ez 3 

where v' = it 2 ij. 

On the other hand, 

1 



dt + -e 3 A -vd z 



A, 



v 



int,0 



[d t v int '°(n, k h ,t, z) + \k h \ 2 v int '°(n, k h ,t, z)] 



inv mt ' (/j, 1 k h ,t,z)e 



+ 



^J-(§-3§)(/*,fcft,t) 




In order that exp(— t/e)L)u™ t + 5u mt '° + v mt '° is an approximate solution 
of (jl.4|) . we project both equations on Ni for I G Z 3 , multiply by exp(iA^), 
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and identify each term. We further apply the following rules in order to 
determine the equations for 5u mt, ° and u™ 1 : 

• all the terms which do not have fast oscillations and are of order 
0(^3^/7) become source terms in the equation on q, 

• all the terms which are either 0(63 or oscillating at a frequency 
1/e become source terms in the equation on 5q. 

We work with a fixed I G Z 3 . Recall that v mt, ° has no purely vertical 
component, i.e. v mt '°(fi, lh, t,z) = if lh = 0. Thanks to formulas (14. 6ft . the 
equation on q reads 
(5.9) 

dtci + \l h \ 2 ci + v'\h\ 2 Cl = -y^jyp [5l{-Xuh,t) - (-l) h 6l(-Xi,l h ,t) 

supplemented with the initial condition 

(5.10) q(0) = (#,|7>, 
and the equation on <5q is 

(5.11) 8^ + (\l h \ 2 + v'\h\ 2 )5c l 



7E 



e ^ 2tt 

i, , n i^M + h*=»\ e i(w>£ 

^°vr|^3 |h| J 

For the time being, we do not specify an initial condition for 5q. Indeed, we 
shall see that it is convenient to choose another condition than —(Ni,v mt, °), 
in order to use the possible decay of S^(fi, lh,t) with respect to t. This choice 
will be made clear in paragraph 15. 41 

As in the previous section, we truncate the large frequencies in <5q. This 
creates an error term in the evolution equation, which is of order 

o- * 1 



eKVy L2 ' 

where K is the truncation parameter, to be chosen later on. We set 

h \h\<K 
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• We now apply to 5u l £ ' the small divisor estimate stated in Lemma [2] 
in the Appendix with 

h\ 3 



s(jj,,l,t) 



%(jM,l h ,t)-{-l) la %(jM,l h ,t) 



d t 5° 3 (Li,i h ,t) - (-i) l *d t 5 l 3 (v,i h ,t) 



7T\l\h 



TT\l\k 



v 5$(ti,l h ,t) - (-l) l ^5l(fi,l h ,t) f f4h\ l h=0 

27T V l3 *°*\l\h \k\ 

from which we deduce that if s < 2, 

ii*4 ; "°«)ii«- < cK i "vsY,{\\4mu+\\4(t)\u} 



\\d s Si(s)\Uds+ sup ||^( 5 )|| 4 
se[o,t] 



, lie intfl II 
+ W du K\t=o\\H s - 



We now choose K such that 



v 1 



K 1 ' 2 



ev, 



i.e. K = e 1 / 2 . We infer that the error term in the evolution equation is of 
order e^V/ 2 in L°°([0, T), L 2 (w)), and that 

intfl i 



\5u 



K IIL-([0,T],^H) < Ce^V/ 2 ^ sup \\8i(t)\\ A 

j te[o,T] 

T 



S JO 



+ 



The operator U is thus defined by 



r intfl || 
\ du K\t=o\\Hi(u)- 



U( T , 3, %)(t) = exp ( --L ) uf\t) + v int '° + 5u™' u , 



intfl 



where u 



int „Antfl X n , int fi 



' are defined by (|5.7p . (j5.6f) and (15. 8p respectively. 



5.2. Choice of the boundary conditions for u BL and u int . We now 
explain how the boundary conditions are chosen. As before, we work with 
kh fixed. Also, since the boundary conditions are all almost-periodic with 
respect to the fast time variable t/e, we work with a fixed frequency fieR. 
Note that this decomposition is allowed by the linearity of the equation. 
We set 

u BL =B(5l5l), 
where the boundary conditions <5°, 5\ are yet to be defined. 
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In order to match the boundary conditions (|1.2p -( fTT3j) with a = 0, we 
must take u BL and u mt such that 

[u BL + ^\ =0 = o{5), 

d z (u BL + ul nt ) lz=1 = o(5), 

K L + O |2=0 = (^), 

(u§ L + u™\ =1 = o(V^S), 

denoting by 5 the order of magnitude of 5°,5 1 , in a sense to be made clear 
later on. 

We now examine each of the boundary conditions independently. 

• At z = 0, the horizontal boundary condition yields 

(5.12) k h ,t) + l,=-x k c k (t) = 0, 

where the vector n(k) is defined in Appendix A (see (16,7p . (16,8p ). 
Since kh is fixed, note that for all fj, £ R, there exists at most one 
&3 G Z such that ^k h ,k 3 = — an d thus the expression above is 
well-defined. 

• Let us now tackle the vertical boundary condition at z = 0. Ac- 
cording to the third section, the vertical component of u BL at z = 
depends on 5®. Precisely, we recall that 



u B L (/i, k h )\ z = = ^fev ^ TTfiikiWxc,! + ik 2 w\o 



(up to exponentially small terms), and 

(a°_,a + ) = p- 1 5° h (fi,k h ). 

As a consequence, in order that the vertical boundary condition at 
z = is approximately satisfied, we choose 

a 

(5.13) (5° = - ^ ^7(^1^,1 + ^2^,2)- 

<T6{-1,1} 

• At z = 1, d z u™ 1 is identically zero by construction of the operator 
U, and thus we infer 8\ = 0. 

• Concerning the vertical component at z = 1, the calculation is the 
same as before. Since (5^ = 0, we deduce that 5\ = 0. 

The above relations (|5.12p - (|5.13p allow us to write 5° in terms of u™ 1 . Con- 
versely, the equation (15. 9p on u 1 ^ depends on 5®, and thus on 5® through 
the operator B. In other words, there is a coupling between the boundary 
condition at the bottom for u BL , and the equation satisfied by u l £ l . Since 
u™* is the only non-vanishing term in L 2 norm, we choose (as is usually 
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done in the rotating fluids literature) to write an explicit equation for u™*, 
and to express u BL in terms of u l £ l . 

5.3. Derivation of the equation for u™ . We now compute the Ekman 
pumping term, that is, the right-hand side in the equation satisfied by c k 
(see (I5.9P ). Notice that if k G Z 3 and ku ^ 0, then |Afc| / 1. In other words, 
the source term in ()5.9[) involves only the part u BL of the boundary layer; 
precisely, with the notations of section O the decay rate of u BL (t, X k , k^) is 

(A±) 2 =i(-A fcT l)+ (l), 

which yields (remember that 3?(A^) > 0) 

7T 



Afc = v 7 ! ± A fc exp (^i-j + o(l). 



Moreover, 



where 



(o°_,a" 



p- 1 [-c fe (t)(n 1 (fe),n 2 (A ; ))] 
-c fe (t)(n_(fc),n + (fc)), 



(n_(fc),n+(A;)) := P^\n x {k), n 2 {k)) 

= ^(ni(k) + in2(k),n 1 (k) - in 2 {k)) + o(l). 

Replacing these expressions in the formula giving #5, we infer 
£5 (-A&, = Cfc(i) ^ ra ^^ (ifciw A T,i + »fc2%,2). 

<J6{-1,1} fc 

We deduce that c& satisfies a linear evolution equation with a damping term, 
namely 

(5.14) ^ + l/c^Cfc + i/|A;3| 2 c fc + ^-A k c k {t) = 0, 

where v 1 = it 2 v and 

^fc — lfc^O^L ^(jfclIi) AJil + *^,2). 

' ' <re{-l,l} fe 

An estimate of 3?(^4fc), where $l(x) denotes the real part of a complex 
number x, is computed in Remark 15.21 below. Using Duhamel's formula, we 
deduce that 

(5.15) |cfc(f)| < exp f-t (\k h \ 2 + u'\k 3 \ 2 + ^U(A k )J\ \{N k ,>y)\. 
We deduce the following Lemma: 
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Lemma 5.1. Assume that 7 £ Vq. Then there exists a unique solution 
v%* G i^ c (R+ , V Q ) n L? 0C ( R + , ^ H) of the equation 



(5.16) 



- A h uf> + [<*] = 0, 



<|t=o = 7. 



where the operator S is defined by 

(5.17) 5 [4 nt ] = ^ A k (N k ,u%*)N k . 

fcez 3 



Hence, in the rest of the section, we take 
(5.18) c fc (t)=7fcexp 



A k t 



k h \ 2 + 



By doing so, we have neglected the vertical viscosity term vd 2 . 

Remark 5.2. (i) Notice that with the scaling we have chosen for the wind- 
stress, there is no Ekman pumping due to the wind. Indeed, the Ekman 
pumping term is of order uj3, which vanishes as e,v — > according to hy- 
pothesis (|4.17p . 

(ii) We emphasize that the operator S constructed above depends on v and 
e through the matrix P, the vectors ui\± and the eigenvalues Xr. However, 
it is useful, for later purposes, to compute the leading order terms in A k , 
which amounts to deriving an equation for the limit of the term u 



■int. 



as e, v 



vanish. Hence we now compute the limit of Ak as e,v — > 0. 

Recall that n\(k) and ri2(k) are given by <\1.8\) . Thus, at first order, 



A k 



\kh\ 
2vr|A;| 



2 E 

<T6{-1,1} 
2 



n\(k) — ian2(k) 



8V2~7r 2 |£f 



1 - At 



= (!-•) + 



= R k + iI k + o{l) 
where Rk and Ik are real numbers given by 



(5.19) 
(5.20) 



R k 



h 



l 



4 



8V2V 2 



1 + A A 



+ 



(ik± 
1 + A fc 



1- A* 



ako 



(l+<) 



8V21T 2 



1 + Afe 1 — Afc 



> 



The Ekman operator appearing in equation (12. 5ft is thus given by the 
following formula, for u €Vq 



(5.21) 



S Ekman [u] ■= ^ ( Rk + U k)( N k, u)N k . 

fcez 3 
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(iii) Recalling the definition of A&, we deduce that 

R k > ci^, 

and thus for every k, for e, v small enough, we have 

R(40 > 

|9f(A fc )| < C. 

To conclude this paragraph, we now give estimates on the boundary con- 
ditions $j in the norm || • 

Lemma 5.3. Assume that 5®, 5® are given by (|5.12p - (|5.13p . Then the fol- 
lowing estimates hold 

\K(t)\\s < (^2\k\ 2{s+1) \c k (t)\^j <C|| 7 ||h. + i, 

and 

\\sim s <\K(t)\\ s+ i<ch\\ H s + 2 {UJ) . 

Proof. The bound on 8® is easily deduced from inequality (|5.15j) together 
with formula (|5,12p and the Cauchy-Schwarz inequality. Concerning the 
other bound, let us recall that if \i = — A/%, for k 6 Z 3 , then the decay rates 
X ± (-X k , k h ) satisfy 

\\±\<c(M + 1 



\kh\ 

Plugging all this estimate into (|5.13p yields the desired inequality. □ 

5.4. Estimates on the boundary layer and corrector terms. Now 

that u™ 1 is rigorously defined by Lemma 15. 11 we may define the other terms 
v mt >°, d u mt fi an d u BL, °. We have gathered in this paragraph some estimates 
which are needed in the proof of Theorem 12.51 

• The boundary layer term of order zero, denoted by u BL , is defined by 

u BL >° = B(SP h ,0), 

where 5® is given by (|5.12p . Thus we deduce that the decay rates A ± (//, kh) in 
the non-resonant part of the boundary layer term u BL, ° are all of order one. 
Consequently, according to (12. 3p . the boundary layer term u BL, ° satisfies 

\\rf Lfi (t)\\ L2{uj) + {ev)-^\\u BL '°(t)\\L^) < C\\6° h (t)\\H H . h )^) 1/4 
(5-22) < C\h\\ mH (e^\ 

Moreover, the definition of c k entails that u BL '° is an approximate solution 
of (jl.4p . with an error term (due to the fact that dtCk does not vanish) 
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bounded in L 2 ([0, T) x u) by 



< 2{ev) 1 ' i ^£ £ (\k h \ 2 + - £ \A k \ 2 ) |^|2 e -2<l^l 2 +Vf^)) dt 

< cM 1/4 (i + ^) 1/2 (^E( 1 + i fc i) 3 i^ 2 ) 1/2 - 



1/2 



The right-hand side of the above inequality vanishes as e, v — > 0, and thus 
the error term satisfies the assumption of Proposition 14.11 

Notice that the Dirichlet boundary condition at z = also generates a 
resonant boundary layer term, namely 

- BL " = ~i E E E 

Ate{-i,i} fc 3 e2Z+i z 3 ez, 3 

sign(Z 3 )/i=l 

We have clearly 

||n BL ' || iO o([ 0iOo)ii 2 ((i ,) ) < C||7|| i 2 (t(;) . 

• The term v int, ° is given by (15. 6|) . in which 5% = and 5§ is defined in 
(|5.13j) . As a consequence, i/"*' satisfies the estimate 

\\v mt '°(t)\\ L i {uj) < C\\5l(t)\\L^ h )^) l/2 
(5.23) < C|| 7 ||^hM 1/2 . 

• At last, the term 5u mt '° is given by equation (|5.1ip . As stated earlier, we 
choose a special solution of (|5.1ip in order to keep track of the exponential 
decay of 5®. Indeed, we have, for all k € Z 3 \ {0}, 

5 3 (-X k ,k hl t) = i%e W ^-(^\k h \ 2 + ^A k ^t^ E ^T-kh-w^. 

Thus we choose for 5cj, |Z| < K, the special solution constructed in Remark 
16.11 in Appendix C. With this choice, we obtain 

|| fox™*' (t) || #2 < 

/ _ \ V2 

Ce^ 2 E 1 il±jM ^ 2l7,| 2 exp (-2^{A k )t) 

Moreover, we recall (see Remark I5.2p that there exists a constant C such 
that |9(Afe)| < C for all k; and in the sequel, we will choose 7 so that 7fc = 
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for k% large enough. In this case, we have 



for e, v small enough and for all k such that % ^ 0. The above estimate 

then becomes 

(5.24) 

_ . 1/2 

5uf fi (t)\\m < C^/V/ 2 ( Y, (! + N) 10 |7 fc | 2 exp L 2 ^(^)t) 

5.5. Conclusion: proof of Theorem 12.51 when a = 0. The idea is to 
use the construction of the previous paragraphs in order to compute an 
approximate solution of the evolution equation (jl.4p . which satisfies the 
boundary conditions up to sufficiently small error terms. We now have to 
quantify the order of approximation required on the boundary condition. 
This is done in Lemma [T] in the Appendix, and thus we build interior and 
boundary layer terms until the conditions of the Lemma Q] are met. 

Let us emphasize that equation (jl.4p supplemented with homogeneous 
boundary conditions at z = and z = 1 is a contraction in I?. As a 
consequence, it is sufficient to prove the Theorem for arbitrarily smooth 
initial data. Thus, without any loss of generality, we assume from now on 
that the initial data 7 only has a finite number of Fourier modes, that is 

7 = ^ kNk ' 

\k h \<N \k 3 \<N' 

Let us now explain the construction in detail. 
• First, we set 

n° :=u int + u BL >°, 

where u int and u BL '° have been defined in the previous paragraphs. We have 
seen that u° is an approximate solution of the evolution equation (II. 4p , with 
error terms which are all o(l) in L 2 . We now evaluate the error on the 
boundary conditions. Indeed, setting 5u := u — u°, we have proved that u 
is an approximate solution of (II. 4p . with some boundary conditions r/°, 77 , 
namely 

Su h \ z=0 = Vh> d z 5u hlz=1 = r/l, 
5u 3 | z=0 = r^, 5u 3 \ z=1 = r)l. 

Thus we have to estimate ^7 := 5u\ t=0 , together with the terms r/°, 77 . 

First, since u™{ t=0 = 7 and u B ^ =Q = 0, we obtain 

(r. o<\ fa - -f, BL >° - v int '° - fa int '° 

where u B ^, w^'q and Sut^^ satisfy the estimates (|5.22p . (|5.23p . and (|5.24[) 
respectively. Thus 

V < C (||7ll^M 1/4 + \h\\H*(ev) l/2 + hWn^v 1 ' 2 ) . 



ROTATING FLUIDS WITH RESONANT SURFACE STRESS 41 

Then, by construction of the operators U and B, the horizontal remainder 
boundary term at z = 1 is exponentially small: indeed, we have <9 z ti^™* =1 = 
0, and consequently, 

(5.26) rjl = -J2 E ""ie^^e^'^e^i 

H,kh ere{-l,l} 

We infer that 

(5.27) \\vlWl < Cexpf-j^) Y, E \S° h (-Xk,k h ,t)\ 2 



< C|ffiexp' ( 



N'y/ev 
Similarly, 



o-e{-i,i} 

and thus 

(5.29) II^Ho < CW'^exp (-^rW) Utt- 

The treatment of the vertical boundary condition at z = is easier. 
Indeed, since S 1 = 0, we have 7/3 = 0, because 

(5.30) V 3 = ~Y1 E aJ^ifch-^e-^e^e^^O. 

There remains to compute r/jj; because of the terms Su^ 1 ' and v mt, °, rf^ is 
the largest term of all. Precisely, we have 

(5-31) rg(t) = -[v^ (t) + 5^ z=0 (t)_ 

i+n \ k h\ 
(5.32) "EE Sc k (t)e ik ^e- iX ^n h (k), 

k h \k 3 \<K 

and thus there exists a constant c > such that for all t > 

Now, the remaining boundary terms 77^, 773, 773 are all of order o(e) accord- 
ing to (|5.27p - (|5.3U|) . Notice furthermore that by construction, 



/ V3 = f° r 3 = 0) 1- 

Juj h 
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Consequently, fh\,n\,n^ all match the conditions of the stopping Lemma HJ 

• We now have to continue the construction with the "bad" part of the 
remaining boundary conditions, i.e. 77?. Let us define the boundary layer 
term 

6u BL >° :=B( V lO). 

According to (|2.3p . 

\\Su BL '°\\ L ~ {[0iOo)>L2H) < Ciasf'Vjo < Ce l ' 2 vV%\\ He , 

and Su BL, ° is an approximate solution of equation (II. 4p with a o(l) error 
term. Moreover, notice that for all t > 0, for all s > 0, 

ll^foWH^K) < C^IMI^exp (-c^) • 
We deduce that for all T > 0, for all s > 

Thus ^u^^fjj satisfies the hypotheses of Lemma [TJ Additionnally, Su< z J® is 
exponentially small, and thus also satisfies the conditions of Lemma [TJ 

• We now define the approximate solution u app by 

u apP ■= u int + u BLfi + 5u BLfi + w, 

where w is defined by Lemma [TJ with the remaining boundary conditions. 
By construction, u app is an approximate solution of the evolution equation 
{T3I), with 

u app\t=0 = U\t=0 + 

and u app satisfies homogeneous boundary conditions at z = and z = 1. By 
a simple energy estimate analogous to that of Proposition 14.14 we deduce 
that 

||^ - u app\\L°°((0,T),L 2 ) ~* VT > 0. 

Since all the terms in u app except u™ 1 and u sm9, ° are o(l) in L? norm, 
Theorem 12.51 is proved. 



Remark 5.4. The proof of Theorem \2.5\ for a = is valid for all ranges 
of e,v such that e, v —* 0. In particular, we do not assume that v = 0(e). 
However, in the case v 3> e, all the modes such that kh 7^ in u 1 ^ are of 
order exp(— cy/u/et), and vanish exponentially for allt > 0. Thus the effect 
of the heterogeneous horizontal modes of the initial data vanishes outside an 
initial layer of size y/e/v. On the other hand, the modes corresponding to 
kh = are not damped, and give rise to resonant boundary layer term u B e ^'° . 
Eventually, for t S> y/e/v, we have 



u (.t) ~ Yl \o,o,k 3 ) N (Ofl,k 3 ) + u 



BL,0 

res 
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6. Towards more realistic models 

6.1. Justification of equation (|1.4p for geophysical models. We now 

explain how our results may give some insight on models of wind-driven 
oceanic circulation, which we recall below. In general, these equations are 
too difficult to deal with in complete mathematical generality, and thus 
crude assumptions are necessary in order to focus on some special phenom- 
ena. Since our aim in this paper is to describe particular kinds of boundary 
layers occurring at the top and the bottom of rotating fluids, we give in this 
regard a few elements on the derivation of the system (jl.4p supplemented 
with (jl.2p - (jl.3p . We emphasize that this derivation is rigourous neither phys- 
ically (since a number of important physical phenomena will be neglected 
in the process), nor mathematically. Our sole purpose is to present some 
motivations for the study of equation (jl.ip . and more generally to derive 
mathematical tools wich may be useful in models of physical oceanography. 

• As a starting point, we recall that the ocean can be considered as 
an incompressible fluid with variable density p; hence, neglecting in a first 
approximation the temperature and salinity variations, the velocity u of the 
oceanic currents satisfies the Navier-Stokes equations, with a Coriolis term 
accounting for the rotation of the Earth 

d t p + u-Vp = 0, 

(6.1) p[d t u + (u- V)w] +Vp = F + puAQ, 

V-u = 0, 

where T denotes as in the first section the frictional force acting on the fluid, 
is the (vertical component of the) Earth rotation vector, and p is the pres- 
sure defined as the Lagrange multiplier associated with the incompressibility 
constraint. 

We assume that the movement to be studied occurs at midlatitudes. At 
such latitudes, we can neglect the variations of the Coriolis parame- 
ter £1 and use the /-plane approximation, which makes the analysis much 
simpler than in the case of the full model. 

The observed persistence over several days of large-scale waves in the 
oceans shows that frictional forces T are weak, almost everywhere, when 
compared with the Coriolis acceleration and the pressure gradient, but large 
when compared with the kinematic viscous dissipation of water. One com- 
mon but not very precise notion is that small-scale motions, which appear 
sporadic or on longer time scales, act to smooth and mix properties on the 
larger scales by processes analogous to molecular, diffusive transports. For 
the present purposes it is only necessary to note that one way to estimate 
the dissipative influence of smaller-scale motions is to retain the same rep- 
resentation of the frictional force 

T = A h A h u + A z d zz u 
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where A z and Ah are respectively the vertical and horizontal turbulent vis- 
cosities, of much larger magnitude than the molecular value, supposedly 
because of the greater efficiency of momentum transport by macroscopic 
chunks of fluid. Notice that A z ^ Ah is therefore natural in geophysical 
framework (see |19j). Moreover, models of oceanic circulation usually as- 
sume that the vertical viscosity A z is not constant (see [31 [18]); we will 
come back on this point later on. 

• Let us now describe the boundary conditions associated with f)6. 1 j) : typi- 
cally, Dirichlet boundary conditions are enforced at the bottom of the ocean 
and on the lateral boundaries of the horizontal domain iVh (the coasts), i.e. 

(g2) u\z=h B (x h )=0 (bottom), 

u\xedu h = (coasts). 

In equation (jl.ip . we have neglected the effects of the lateral boundary 
conditions by considering the case when Uh is the two-dimensional torus. 
Of course such an assumption is not physically relevant. It is well known for 
instance that the lateral boundary layers, called Munk layers, play a crucial 
role in the oceanic circulation, in particular in the western intensification of 
currents. Moreover, for the sake of simplicity, we did not take into account 
the topography of the bottom in (jl.3p . The topographic effects described 
by the function Kb should actually modify the Ekman boundary layer and 
consequently the limit equations, even if the variations of the bottom are 
small (see [5] and [8] for instance). 

We assume that the upper surface, which we denote by T s , has an equation 
of the type z = hs(t,Xh)- As boundary conditions on T s , we enforce (see 

®) 

-Q^-0<z<h s {t,x) + ^ v x{i-0<z<h s (t,x h ) u ) = U 

where £ is the total stress tensor of the fluid, and a w is a given stress tensor 
describing the wind on the surface of the ocean. In general, T s is a free 
surface, and a moving interface between air and water, which has its own 
self consistent motion. In (11.21). we have assumed that 



h s (t,x h ) = D, 

where D is the typical depth of the ocean. Hence (|1.2|) is a rigid lid approx- 
imation, which is a drastic, but standard simplification. The justification of 
(|1.2j) starting from a free surface is mainly open from a mathematical point 
of view; we refer to PQ for the derivation of Navier-type wall laws for the 
Laplace equation, under general assumptions on the interface, and to [13] for 
some elements of justification in the case of the great lake equations. Nev- 
ertheless, from a physical point of view, the simplification does not appear 
so dramatic, since in any case the free surface is so turbulent with waves 
and foam, that only modelization is tractable and meaningful. Condition 
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(|1.2p is a simple modelization which already catches most of the physical 
phenomena (see [19]). 

• Let us now evaluate the order of magnitude of the different parameters 
occuring in (16. ip . and write the equations in a nondimensionalized form. 
First, since the variations of density are of order 10 -3 , we neglect the effects 
of the variations of p in (|6.ip and we assume that 

p = po = 10 3 kg • m~ 3 . 

Moreover, we set 

u h = Uu' h , u 3 = Wu' 3 , 

Xh = Hx' h , z = Dz' , 

where U (resp. W) is the typical value of the horizontal (resp. vertical) 
velocity, H is the horizontal length scale, and D the depth of the ocean. In 
order that u'(x') remains divergence- free, we choose 

UD 

W = . 

H 

Typical values for the mesoscale eddies that have been observed in western 
Atlantic (see for instance [19J) are 

U ~ 1 cm • s _1 , H ~ 100 km, and D ~ 4 km. 

With these values, we get 

- 10 -3 , 



HQ 

and hence t«l (notice that the parameter e is dimensionless) . Thus the 
asymptotic of fast rotation (small Rossby number) is valid. 

A typical value of the horizontal turbulent velocity is ~ 10 6 kg • m" 1 • s~ 
(see [3]), which yields 

Ah 

~ i. 

PoUH 

In general, the vertical eddy viscosity A z is not assumed to be constant; in 
[31 [18] , the authors consider a vertical viscosity which takes the form 

a ( , (a 5 9d z p ' 

A z = po \v b + 1/ 1 



Po\d z Uh\ 



and they assume in their numerical computations that A z > 1 kg • m 1 • s 1 . 
The quantify 

(6.4) Ri := -(gd z p)/( Po \d zUh \ 2 ) 

is called the local Richardson number. Equation (jl.4p corresponds to a con- 
stant approximation for the viscosity A z ; this is largely inaccurate, since 
according to [18], measurements show that the value of A z is usually large 
inside the boundary layer (say, 3 to 10 kg • m" 1 -s -1 ), but substantially 
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smaller in the interior (under the thermocline). However, since we are pri- 
marily interested in the boundary layer behaviour, we only retain the typical 
boundary layer value A z ~ 5 kg • m -1 • s^ 1 , which yields 

HA, 



~ 5 • 10" 3 . 



• poUD* 

Hence we also have v -C 1, which justifies our assumption of vanishing 
vertical viscosity. Notice that the parameter v is also dimensionless. 

Thus the nondimensionalized system (see for instance jl9j, ITO] ) becomes 

, v d t u' + v! ■ Vu' + ~e 3 A v! + ( 7o P / J - A fc u' - vd zz v! = 0, 
(6.5) e \jiOzP J 

V • v! = 0, 

where 5 := D/H is the aspect ratio. The boundary conditions are (jl.2|) . 
(fL3l) . with 

The equation for the boundary layers at z = 1 and z = in the above 
system is exactly the same as in (|1.4p . Thus, we believe that the phenomena 
we have highlighted (atypical size of boundary layers for resonant forcing, 
possible destabilization of the fluid for large times) may prove to be useful 
when studying models of oceanic circulation. However, we do not claim 
that our results truly apply as such to realistic geophysical models, since, 
as mentionned above, a series of drastic simplifications have been made. 
Furthermore, some assumptions of Theorem 12.51 such as ()4.17j) . are purely 
technical, and do not have any physical ground. Thus, we now turn to some 
possible mathematical extensions of Theorem 12.51 to more realistic models. 

6.2. Possible extensions. The previous study allows to characterize the 
linear response of a rotating incompressible fluid to some surface stress, 
which admits fast time oscillations and may be resonant with the Coriolis 
force. In addition to the usual Ekman layer, we have exhibited another - 
much larger - boundary layer, and a resonant boundary layer term, the size 
of which depends on time. Note that these effects do not modify the mean 
motion (i.e. the 1? asymptotics) when considering moderate times, say for 
instance t <C — . 

• Extensions to nonlinear equations. In order to take into account 
more physics in our model, the first point is to understand the nonlinear 
response of the fluid to the same surface stress. In other words, we are in- 
terested in the asymptotic behaviour of the full Navier-Stokes-Coriolis equa- 
tion (|6.5p - (ll.2p - (|1.3p including in particular the nonlinear contribution of the 
convection. 

In the case of a non-resonant forcing, the asymptotic motion of the fluid is 
obtained by some filtering method : there is indeed two time scales, a rapid 
time scale at which the fluid oscillates according to the modes of the linear 
penalization, and a slow one which characterizes the nonlinear evolution 
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of the wave enveloppes. The boundary effects do not play any role in the 
nonlinear process since they are localized in the vicinity of the surface. They 
contribute to the envelope equations only by the Ekman pumping. In the 
case of a resonant forcing, the boundary effects - which are not expected to 
be localized in the same way - could play a different role. 

• Towards more physically relevant models. The present theory of 
the wind-driven circulation of a fluid of uniform density is actually inade- 
quate to capture the velocity structure of the oceans. We indeed expect the 
wind forcing to modify in depth the circulation. The profile arising from the 
resonant part of the forcing and the Ekman pumping are not enough to get 
a relevant description of that vertical structure. 

We will mention here many phenomena that have been neglected in our 
study and which seem to be crucial to obtain realistic models. 

(i) we first need to consider the variations of the Coriolis parameter, 

keeping at least the /3-plane approximation : 

n = f + p y 

where y is the coordinate measuring the latitude. Such a spatial de- 
pendence of O is necessary to derive Sverdrup's theory of horizontal 
transport, which is still one of the foundations of all theories of the 
ocean circulation (see [20] for instance). 

From a mathematical point of view, we refer to [5] [7] and references 
therein for some preliminary studies on inhomogeneous rotating fluids. 

(ii) the vertical structure of the ocean cicrulation is also related to the vari- 
ations of the density p, the so-called stratification of the oceans. The 
theoretical works of Rhines and Young [21] have brought some under- 
standing about geostrophic contours, potential vorticity homogeneiza- 
tion and their role in shaping the pattern of circulation. Luyten, 
Pedlosky and Stommel [15] have then developped a theory for the 
full density and velocity structure of the wind-driven circulation by 
going beyond the quasi-geostrophic approximation to consider the im- 
portant effect of the ventilation of the thermocline which occurs as 
oceanic density surfaces rise to intersect the oceanic mixed layer. 

However, to our knowledge, there is no mathematical contribution 
on that topics, the first difficulty being to determine some suitable 
functional framework to deal with the inhomogeneous incompressible 
Navier-Stokes equations. Moreover, the behaviour of the fluid is ex- 
pected to depend in a crucial way on the order of magnitude of the 
Richardson number Ri, defined in (16, 4j) above: when Ri is small (say, 
Ri < 1/4), instabilities may develop, leading in turn to some turbu- 
lent mixing across layers of equal density. We refer to [24] for more 
details. 

(iii) we finally have to take into account the bottom topography which may 
have an important contribution to the mean circulation as proved for 
instance in [5] or [HJ. 
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The crucial point to understand these features from a mathematical point 
of view is to get a description of the boundary layer operator which is not 
based on the Fourier transform, but on the spectral decomposition of the 
Coriolis operator. The Coriolis penalization becomes indeed in the two first 
cases a skew-symmetric operator with non-constant coefficients (depending 
on f2 and p). We therefore have to develop new tools to obtain the asymp- 
totic expansions in a more abstract and systematic way. 

Appendix A: spectral results on the Coriolis operator 

For the sake of completeness, we recall here - essentially without proof - 
some fundamental properties of the Coriolis operator leading to (|1.8p . For 
a detailed study of these spectral properties we refer for instance to [3]. 

Extending any u € Vq on [— 1, 1] x T 2 as follows 



(6.6) 



Uh(xh, z) = Uh(xh, -z) and u 3 (x h , z) = -u 3 (x h , -z) 



(which is compatible with the incompressibility constraint V • u = 0) we 
obtain a periodic function, so that it is possible to use some Fourier decom- 
position. 



Setting 



(6.7) 



and 



(6i 



n\{k) 
n 2 (k) 
n 3 (k) 



1 



2ir\k hi 
1 

2n\k h \ ' 



(ik 2 + fciAfe) 
-iki + k 2 \k) if k h / 0, 



27VN 2 + (^ 3 ) 2 
sign(/c 3 ) 



m(k) 

n 2 {k) 
n 3 (k) 



2tt 



else, 



2tt 




what can be proved actually is that the family (N^) defined by 
N k = exp(ik h ■ x h ) 



ni(k) cos(irk 3 z) 
n 2 (k) cos(-Kk 3 z) 
n 3 (k) sin(irk 3 z) 



is an hilbertian basis of Vq constituted of eigenvectors of the linear penal- 
ization, satisfying (|1.8|) . 

Appendix B: the stopping condition 

We have postponed here the statement and the proof of the stopping 
condition since it is just a technical result (based on straightforward com- 
putations) which is used in several places (Sections 4 and 5). 
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Lemma 1 (Stopping condition). Let S ^ 1 G L°°(R + , H 3 (ujh)) be two fam- 
ilies such that 



J (Si - S° 3 )dx h 





and 

-W\\w(u h ) ~> 0> ll**llHS( Wh ) -»• and \\d t S l \\ H i M -> as e -> 
XTien i/iere exists a family w G L°°(R + , L 2 (S7)) uref/i V • u; = stic/i t/iai 

u>|*=o = 5°, w 3 \ z=1 = S\ and d z w h \ z=1 = 5\ 
and satisfying the following estimates 



MIl 2 (h) 



and 



dtw H — Lw — vd zz w — A^ui 
e 



— > as e — ► 0. 



Proof. Here we have to build a family «; G L°°(R + , L 2 (S7)) with V ■ to = 
such that 

w\ z=0 = 5°, w 3 \ z=1 = Si and d z w h \ z=1 = S\. 

Of course it is not uniquely defined. We just want to obtain one such family 
satisfying further suitable continuity estimates. 

Given any two-dimensional vector field Wh, we get a divergence- free vector 
field by setting 

w 3 (x h ,z) = w 3 (x h ,0) - / (dxwi + d 2 w 2 )(x h ,z')dz' . 

Jo 

In order that the boundary conditions on w 3 are satisfied, the only condition 
on Wh is therefore 

f\d lWl + d 2 w 2 )(x h , z')dz' + 5\{x h ) - S° 3 (x h ) = . 

JO 

We therefore choose 

wi(x h , z) = 5i(x h ) + S\(x h )z + di<f>(x h )z(l - z) 2 , 
w 2 (x h ,z) = 5%(x h ) + S\(x h )z + d 2 4>(x h )z(l - z) 2 , 

with 

V ft • 5° h + l -V h ■ 5\ + ^A^ + Sl -5° 3 =0. 

Standard elliptic estimates give for any s > 

WHh^m < C(\\^\\ H . (uh) + II^H^K))- 

Therefore 

\H\h*(Q) < C(\\S°\\ H 3 {oJh) + H^Htfs^)) 
so that, using the assumptions on 6°, S 1 , 

IMItf 2 (fi) ^ as e ^ 0. 
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Furthermore, since w is given in terms of J ,^ 1 by linear relations with 
constant coefficients, 

\\dtw\\ L2{Q) < C(\\d t 5°\\ H i M + II^Hhi^))- 

We conclude, using again the assumptions on S ^ 1 that 



dtw H — Lw — vd zz w — AhW 



as e -> 0. 



L 2 (n) 



□ 



Appendix C: the small divisor estimate 

We recall here the by-now standard arguments used to obtain some es- 
timate for the solution to fast-oscillating linear equation with non-resonant 
source terms : 

(6.9) d t w + -F(w) - uA h w - vd 2 zz w = E 

where the horizontal Fourier mode lh is fixed and 

ft k 3 ez 

We further assume that the frequencies \i belong either to {—\i,h G Z 3 }, 
or to some finite set M. 

The small divisor estimate is the following: 

Lemma 2. Let w be the solution of A6.9\) . i.e. for all I = (lh, I3) with I3 £ Z, 

d tWl + {\l h \ 2 + v'\h\ 2 ) Wl = Y, s(iu,l,t)e^ +x '^ 

where v' = it 2 v. 

Then there exists a constant C such that for all t > 0, r > 0, for all 
K > 0, we have 

\\F K w(t)\\ H r (lj) < Ce{\\s ]t=0 \\ r , K e^p(-{\l h \ 2 + u'l 2 )t) + \\s(t)\\ r , K } 
+ Ce [ \\d u s{u)\\ 7 , K ^ V {-(\l h \ 2 + u'l 2 )(t-u)) du 



+ Ce sup || s(it) I | ri _K- 
ue[o,t] 

+ WKW\ t =o\\H r {u))i 

where the norm \\ ■ \\ r k is defined by 



Ht)\\r,K ■■= E E |*»| 8 |iri*(-Afc,i > t)| 5 



\1\<K fc 36 z 



+ E E (i + iH=iRI 4 )RrK^U)l 2 - 
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We recall that the notation ¥k stands for the projection onto the vector 
space generated by Nk for \k\ < K. 

Proof. For all K > 0, define 

w K ■= ¥ K w = WlNi. 

\k\<K 

We deduce from DuhameFs formula that 
\ Wl (t)\ < \ Wl (0)\ exp(-(\l h \ 2 + v'll)t) 

(6.10) + f V s(/i,/,n)e^ +Ai )texp(-(|^| 2 + i/ , /|)(t-u)) ( iu 

Integrating by parts, we get 

s(/i, I, u)e i{ - Xl+ ti-* exp(-(|Z fc | 2 + u'l 2 )(t - it)) du 
Z, t)| + tt-^— r|s(/x, Z h , 0)|e-^l 2+ ^)* 



< 



+ 



+ 



\h+v\ 
e 

|A*+/u| 

c 



|Ai + mI 

KKfcl 2 + v'\k\ 2 )\s(v,l,u)\ exp(-(|Z h | 2 + y'll){t - u)) du 



\d u s(n,l,u)\ exp(-(|Zh| 2 + v'll)(t - u)) du. 



|Aj + mI 

Plugging this inequality back into (|6.10p . we deduce that 

\wi(t)\ < \wi(0)\exp(-(\l h \ 2 + i/ll)t) 

\s(fi,l,t)\ 



+ Ce y: 



|Aj + mI 



+Ce ^ l^^yi e xp(-(|Z fe | 2 + ^Z|)t) 



rt 



\k+fi\ 



+Ce I F^expi-Qlh]* + v%)(t-u))du 



where 



Fi(u) 



£ ivb^ uS(/i ' u)l 



Mt^-A; 



+ (N 2 W|z 3 | 2 ) £ iy^I^U)!- 

There remains to derive bounds for quantities of the type 
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Remember that either [i = — A& for some k = (l h ,h) e Z 3 with k 3 / -Z 3 , 
or /i G M, where M is a finite set. Thus 



(j : ^ K "' ^ •'' ,l ) 
£ 2 fe.ra |s( - A " v "^) 

+2 (s^ KMu)l ) 

^ C E l fc 3| 2 |x \ . 2 |g(-A fc ,f,n)| 
' — \M-*k\ 



Notice that the function £3 1— > A/ is monotonous for 1^ fixed. Hence |A/ — A^j 
is minimal for k 3 = Z3 ± 1. Consequently, is is easily checked that for all 
h € Z, 

ml 

Moreover, if fi G M, then either /x ^ {0,1,— 1}, and in this case 

1^ -^r 1 <C, 

or /i = 0, and then 

|A, - Ml" 1 < C^L, 

or /x G {1, —1}, and then 



N 2 ' 

Gathering all these results we get 
\w t (t)\ < \ Wl (0)\+CeD°(t) 

(6.11) +Ce J D h u ) ex p(^- + du ' 
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where 



Y,\ k 3\ 8 \ s (- x kJ, 



1/2 



exp(-(\l h \ 2 + u'll)t) 



+ Yl (l + lH=l|Z| 2 )KM,^0)|e X p(-(^| 2 + ^|)t) 

1/2 



Jl^-A; 



+ 



^|A:3| 8 | S (-A fc ,/,t)p 



fc 3 



+ ^ (i + i H=1 |/| 2 )|^,U)l 



and 



^|A: 3 | 8 |^(-Afc,/,n)| ; 



1/2 



+ E E (i + iH=iif)i^(/i,/,n)i 



M^-A, 



+ (|/h| 2 +^K3 



El^| 8 |s(-A fe) Z,«)p 



1/2 



+ (|Z h | 2 + i/|i 3 | 2 ) E ( 1 + 1 M=il i l 2 )l s 4 '^«)l- 



ht^-a. 



The estimate of Lemma [2] follows. □ 

Remark 6.1. Assume that the Fourier coefficients of s have exponential 
decay, meaning that for all (fJ-,1), there exists sq([jl,1) £ C, and c(fi,l) G C 
with nonnegative real part such that 

s(fjt, I, t) = soi/x, I) exp(-c(/i, l)t). 

Then provided the sequence so(fM,l) is sufficiently convergent, a special so- 
lution of (|6.9p can be built, which preserves the exponential decay property. 
Indeed, for all I £ Z 3 , set 

exp (i(Xi +n)l - c(fj,,l)t) 



c(n,l) + \l h \ 2 + v\h\ 



Then it can be readily checked that w is a solution of (16. 9p . and moreover 

1 



K*)l < e E 



|A, + M-e9f(c(M,0)l 



I «o (a*, 01 ex P (-^( c (a*> 0)*) ■ 
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